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Abstract 



We study a nonlinear ground state of the Gross-Pitaevskii equation with a parabolic 
potential in the hydrodynamics limit often referred to as the Thomas-Fermi approximation. 
Existence of the energy minimizer has been known in literature for some time but it was only 
recently when the Thomas-Fermi approximation was rigorously justified. The spectrum of 
linearization of the Gross-Pitaevskii equation at the ground state consists of an unbounded 
sequence of positive eigenvalues. We analyze convergence of eigenvalues in the hydrodynamics 
limit. Convergence in norm of the resolvent operator is proved and the convergence rate is 
estimated. We also study asymptotic and numerical approximations of eigenfunctions and 
eigenvalues using Airy functions. 

1 Introduction 

Recent experiments in Bose-Einstein condensation has stimulated an intense research around 
the Gross-Pitaevskii equation with a parabolic potential [PS]. Considered in a one-dimensional 
cigar-shaped geometry and in the limit of a compact Thomas-Fermi cloud, the repulsive Bose gas 
is described by the Gross-Pitaevskii equation in the form 

iut + £ 2 u xx + (1 — x 2 )u — \u\ 2 u = 0, (1-1) 

where u = u(x,t) is a complex- valued amplitude, the subscripts denote partial differentiations, e 
is a small parameter, and all other parameters are normalized to unity. 

Existence of the ground state u = %(x) for a fixed, sufficiently small e > 0, where r] e is a 
real-valued, positive-definite, global minimizer of the Gross-Pitaevskii energy 

E e {u) = j Qe 2 K! 2 + \{x 2 ~ 1)M 2 + \\u\ A ^j dx 

in the energy space 

Hi = {u e H 1 ^) : xueL 2 (R)}, 

has been proved in the literature long ago (see, i.e., Brezis &: Oswald [BO] ). Recent works of 
Ignat & Millot [IM] and Aftalion, Alama, & Bronsard [AAB] have focused, among other problems 
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related to existence of vortices in a two-dimensional rotating Bose-Einstein condensate, on the 
rigorous justification of the Thomas-Fermi asymptotic formula 

, , f (l-x 2 ) 1/2 , for M < 1, n 
= 1 0, for \x\ > 1, (L2) 

which was believed to be a weak limit of r/ e (x) as e — > since the work of Thomas [T] and Fermi 
[F] . To be precise, Proposition 2.1 of [IM] and Proposition 1 in [AABJ state that r\ e {x) converges 
to t]q(x) as e — > in the sense that 



(1 " Ce 1 /*) < < 1 for \x\ ^ 1 - e 2 / 3 

< r/ £ (x) < Ce 1 / 3 exp for \x\ > 1 - e 2 / 3 , 



(1.3) 



for an e- independent constant C > 0. (The results of |IM1 lAABj are formulated in the space 
of two dimensions, but the extension to the one-dimensional case is trivial.) It was proved in 
|IM] that \\rj e — f/ollc 1 ^) ^ Cr£ 2 for any compact subset K C (—1, 1), which justified the WKB 
approximation of the ground state considered earlier by formal expansions (see, i.e., [BK]). 

We are concerned here with the spectrum of linearization of the Gross-Pitaevskii equation 
(jl.ip at the ground state rj e , which is defined by the eigenvalue problem 

- e 2 u" + (x 2 - 1 + 3i] 2 )u = -Xw, -e 2 w" + (x 2 - 1 + i] 2 )w = Xu, (1.4) 

where (u+iw)e xt + (u— iw)e xt is a perturbation to %. The eigenvalue problem (jl.4|) determines the 
spectral stability of the ground state rj E with respect to the time evolution of the Gross-Pitaevskii 
equation (jl.ip and gives preliminary information for nonlinear analysis of orbital stability and 
long-time dynamics of ground states. More complex phenomena of pinned vortices (dark solitons) 
on the top of the ground state can also be understood from the analysis of eigenvalues of the 
spectral problem (jl.4jl (see, i.e., |PK] ) . 

In what follows, we shall simplify the spectral problem (jl.4p and replace r] e by rjQ. We do not 
claim that eigenvalues of these two problems are close to each other but, given a complexity of 
the problem, we would like to deal with a simpler problem in this article. Therefore, we analyze 
here solutions of the model eigenvalue problem defined explicitly by 

— e 2 u" + 2(1 — x 2 )u = — Xw, —e 2 w" = Xu for \x\ < 1, - . 

-e 2 u" + (x 2 - l)u = -Xw, -e 2 w" + (x 2 - l)w = Xu for \x\ > 1, ^ ' ' 

with appropriate matching conditions at x = ±1. It will be left for the forthcoming work to 
study solutions of the original eigenvalue problem (jl.4p with rj e = r/o + O^oo^^e 1 / 3 ), according 
to the bound (11. 3p above. 

Formal weak solutions of (|1.5p have been constructed in the pioneer work of Stringari [S] 
and have been used in a more complex context of three-dimensional anisotropic repulsive Bose 
gas in [FCSG1 IEGQ| . To recover these solutions, let us denote A = ie^ 1 / 2 and drop —e 2 u" 
term in the first equation of (|1.5p . Then, the model eigenvalue problem is closed at the singular 
Sturm-Liouville problem 

- 2(1 - x 2 )w" = jw, -Kx<l, (1.6) 

which has a C 2 solution on [—1, 1] for 7 7^ if and only if w(l) = w(— 1) = 0. We will show in 
Lemma 13.41 below that the only solutions of (|1.6p with w(l) = w(— 1) = are the Gegenbauer 
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— 1/2 

polynomials w(x) = C n+1 (x), which correspond to eigenvalues at 7 = j n = 2n(n + 1), where 

n ^ 1 is an integer. Solutions w(x) = C n+l (x) of (jl.6p on the interior domain [—1,1] are 
completed with the zero function w = on the exterior domain |x| ^ 1. In this way, we glue 
together weak solutions of system f)1.5|) in the hydrodynamics limit e = 0. It is the main goal 
of this article to develop a rigorous justification of persistence of eigenvalues {7 n }neN for small 
non-zero values of e. Our main result is the following theorem. 

Main Theorem. Spectral problem U.5\) for e > has a purely discrete spectrum that consists 
of eigenvalues at A = ±ie(7 n e ) 1 / 2 ; where the set {7n, e }neN * s sorted in the increasing order 

while 

7n, £ — > In as e ^ 
for every fixed n £ N. Moreover, for any fixed 5 > 0, i/iere exists C n > sttc/i i/iai 

|7n, e -7n| 

/or sufficiently small e > 0. 

Remark. 27ie convergence rate of eigenvalues is not sharp and our numerical results indicate 
that the convergence rate is 0(e 2 ) for a fixed n G N. 

Before going into technical details of our analysis, we mention three relevant applications 
where eigenvalues of the singular Sturm-Liouville problem (|1.6p have appeared recently. 

• Propagation of self-similar pulses in an amplifying optical medium is described by the 
Gross-Pitaevskii equation with a parabolic potential [BTNNJ 

iU T + t- 2 U^ + (1 - i 2 )U - \U\ 2 U = 0. 

The small parameter e = t~ 1 changes with the time r due to evolution of the self-similar 
optical pulse in the presence of the gain. The decomposition of perturbation to the optical 
pulse via Gegenbauer polynomials is used for understanding the effects of higher-order 
dispersion and gain terms on the long-term optical pulse dynamics [BT| . 

• Analysis of radiation from a dark soliton oscillating in a wide parabolic potential was studied 
in |PFKj using asymptotic multi-scale expansion methods. The analysis leaded to the wave 
equation with a space-dependent speed 



Urr={(l-eW 



Eigenvalues of the wave equation are given by eigenvalues of the Sturm-Liouville problem 
(jl,6p . The corresponding eigenfunctions are needed to match the dark soliton with its 
far-field radiation tail and to predict radiative corrections to the soliton dynamics [PFK]. 

Numerical approximations of eigenvalues of the spectral problem associated with a dark 
soliton in the Gross-Pitaevskii equation 

iU T + % + i 2 )U - \U\ 2 U = 
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showed convergence of eigenvalues in the limit \x — > oo [PK] . It was observed that the whole 
spectrum consisted of eigenvalues associated with the ground state and an additional pair 
of pure imaginary eigenvalues. The countable infinite set of eigenvalues associated with 
the ground state corresponds to the set of eigenvalues of the Sturm-Liouville problem (|1.6j) 
after an appropriate rescaling transformation of £, r, and U. 

This article is organized as follows. Section 2 discusses properties of the two Schrddinger 
operators that define the spectral problem (|1.5|) as well as the properties of their product. Sec- 
tion 3 gives a proof of the Main Theorem. Section 4 is devoted to asymptotic and numerical 
approximations of eigenvalues of the spectral problem (II. 5h . In the Appendix, we give the proofs 
of several technical lemmas used in the article, as well as the description of the numerical method. 

Notations. In what follows, if A and B are two quantities depending on a parameter p in a 
set V, the notation A(p) < B(p) indicates that there exists a positive constant C such that 

A(p) ^ CB(p) for every p E V. 

The notation A(p) ps B(p) means that A(p) < B(p) and A(p) > B{p). We say that a property is 
satisfied for < e <C 1 if there exists £o £ (0, 1) such that the property is true for every e E (0,£o). 
If E and F are two Banach spaces, C(E, F) denotes the space of bounded linear operators from 
E into F, endowed with its natural norm 

\\ u \\c(E,F) = SUp — : . 

If E = F, we simply denote C{E) = C(E, E). The dual space of E is denoted by E' = £{E, R). 
If S is a subset of R, I5 denotes the characteristic function of S: 



l if xeS, 
if x <£ S. 



If / is a function defined on some set D and S C D, f\g denotes the restriction of / to the set S. 
Finally, B L 2 denotes the unit ball of L 2 (R). 

2 Preliminaries 

2.1 The operator U_ and its inverse 

Let L e _ be the Friedrichs extension of — d% + p e {x) on L 2 (R) for e > and 

p e {x) = ^(x 2 - 1)1 {N>1} . 

Since p e (x) ^ for any x £ R, L £ _ is a positive self-adjoint operator. Since p £ {x) — > +cxd as 
x — > 00, LI has compact resolvent. The domain of Li, 

D(L £ _) = {fe L 2 (R) : -d 2 x y + p £ v E L 2 (R)} = {v? E # 2 (R) : x 2 ^ E L 2 (R)} =: H 2 , 
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is contained in its form domain 

Q(L £ _) = {(/?£ ff X (K) : xp G L 2 (M)}. 

If y € L>(Li) is in the kernel of Li, then J K (|<9 X (/9| 2 + p e |(/?| 2 ) dx = 0, which implies <p = 0. 
Therefore c(Li) and Li is invertible. In the following lemma, we state that the inverse of 
Li is uniformly bounded in C(L 2 ) as e — > 0. 



Lemma 2.1 Lor 0<e<l, 



\\(Ll 



|£(L3) 



1. 



Proof. See Appendix IA.ll ■ 
Using Lemma l2,lt we give estimates on various norms of (Li) -1 for sufficiently small e > 0. 
Lemma 2.2 For < e -C 1, 



9 X (L1) -1 £(L 2 (R)) 


< 


1, 


(2.1) 


ll 1 {|o;|>l}^(- t '-) _ £(L 2 (R)) 


< 


1/3 
c ) 


(2.2) 


l 1 {|a;|>l}(- t 'i) - l/;(L2(R)) 


< 




(2.3) 


|| ^ (LI) 1 \\c(L 2 (R),L°°(R)) 


< 


1, 


(2.4) 


1{|cc|>1}(-^-) 1 \\c(L 2 (R),L°°(R)) 


< 


e 2 l\ 


(2.5) 



Proof. Let us take e > sufficiently small, / 6 L^2, and denote 93 = (Li) x /. By Lemma [2. 1| 

IIvI|l 2 (r) < 1- (2-6) 

Moreover, ip satisfies the second-order differential equation 

-p"+p £ p = f, x£R. (2.7) 

Multiplying (|2.7p by 93, integrating over R, using the Cauchy-Schwarz inequality and (|2.6p . we 
get 



/ \<p'\ 2 dx+ / p e |v?| 2 dx= / ftpdx < ||/Hl2(r)||vIIl2(r) ^ 1, 

JR 7|a|>l JR 



(2. 



which directly proves (|2.ip . Proceeding like for ()2.8|) . but integrating on [l,+oo) instead of M, 
we obtain 



+00 



+00 



</?'| dx+ / p £ \p\ dx ^\ip(l)\\ip'(l)\ + \\<p\\ L 2 {lj+00 y 



(2.9) 



Then, we observe 



MI L 2( 1+£ 2/3 )+00 ) 



+OO 



l+ £ 2/3 X 2 - 1 



p £ \tp\ 2 dx 



(l + e 2 / 3 ) 2 -Ui +£ 2 



2/3 



Pe|<£>| 2 cte 



< p 4 / 3 



(2.10) 



Since ip" = —f on (—1,1) and thanks to bound (|2.ip . Sobolev's embedding of H (—1,1) into 
L°°(-l,l) yields 

II^IU«.(-i,i) < y\\m(-i,i) £ Wh*(-i,i) + II/IIl^-i,!) £ i- (2-ii) 

The triangle inequality yields 

(1,1+e 2 / 3 )- (2.12) 

By the Taylor formula and the Cauchy-Schwarz inequality, 

,+oo)- 

(2.13) 

Let us introduce the new variable £ = (x — l)/e 2 / 3 and the function <p(£) = tp(l + e 2 / 3 £). Then, 

\m\m ( i, + oo) = ^ /3 ii^iii 2(1+£2 /3 i+00) + ^ 2/3 ii^ni 2(1+£2 /3 )+00) (2.i4) 

Thus, by Sobolev's embedding of -ff 1 (l,+oo) into L°°(l,+oo), (|2.14l) provides the bound 

\<p(l + e 2 / 3 )| = |^(1)| < e 1/3 ||^ll i2(1+e2 /3, +oo) + e" 1/3 |M| L 2 (l+e 2 / 3 ,+oo)- (2.15) 
Concatenating dZTDjl . QZty . (|23I]) . (|2TT2|) . (jmjl and <|27r5|> . we obtain 

ll^lli 2 (l, + oo) + ^l^lli 2 (l +£2 /3, +oo) < ^ /3 |l^llL 2 (l, + oo)+£- 1/3 |^ll i2 (l +e2 /3, +00 ).(2.16) 

There exists C > such that (|2.16p can be rewritten in the form 

(ik'lUwoo) - Ce 1 / 3 ) 2 + ^ (lMLw/3, + oo) - Ce) 2 < £ 2/3 . 
Therefore, ||</?'||l2(i,+oc) < e 1/3 and |M| L 2 (1+e 2/3 ]+00 ) < e. Using also (|2.13|) and (|2.15|), we deduce 

II^IL2(i,i+ £ 2/3) < e 1/3 ||vlL<»(i ) i+ e a/8) ^ e, 

and thus |M| i2 ( lj+00 ) < e. Similar computations on (— oo, — 1] complete the proof of (|2.2h and 
Sobolev's embedding of H 1 ^) int o L°°(R + ) for = 99(1 +e 2 / 3 £) yields 

l|y]U°° (l,+oo) = [|<PlU°°(K+) ^ IMItf 1 (K + ) ^ II^'IIl 2 (R + ) + II^IIl 2 (r + ) 

< e 1/3 |l^llL 2 (i, + oo) + ^ 1/3 ||^llL 2 (i, + oo) <e 2/3 . (2.17) 

Combined with a similar estimate for ||y||i<x>(_oo,-i)) we get (|2.5p . Finally, Sobolev's embedding 
of into L°°(M + ) for = e 2 / 3 yj'(l + e 2 / 3 £) similarly yields 

II^IIl°°(i,+oo) < ^ 1/3 II^"IIl2 (1>+oo) +e~ 1/3 ||v? / || i2(li+oo) . 

Therefore, the bound (|2.4p holds if Hv?" Hz, 2 (1,00) ^5 e -1 ^ 3 since ||v?'llz,°°(-oo,-i) is estimated similarly 

and || V 7 ' || z,°° ( 1,1) ^ s gi ven by the bound (|2.11|) . Since if £ D(L e _) = H2, lim a; _*ooPe¥ , ¥' / = 0, and 

the bound ||¥>"||z, 2 (i,oo) ^5 e~ 1 ^ follows from integration by parts: 

/+00 r+00 r+00 

Wfdx - 2 J VeW "dx + J p 2 ip 2 dx 

/+00 r+00 r+00 

(<p") 2 dx + 2 J p e (ip') 2 dx + J p 2 cp 2 dx 

2 r+00 2 
--/ ^ 2 dx--<p 2 (l), (2.18) 
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where the second and third terms in the right-hand-side are positive and the last two terms are 
estimated from (I2.3[) and (12. 5|) . ■ 



2.2 The operator L £ + and its inverse 

Let L £ + be defined similarly to Li as the Friedrichs extension of —d 2 + q e (x) on L 2 (IR) for e > 0, 
where ^ 

?e(aO = ^2 [ 2 (* - a;2 ) 1 {|a:|<l} + ~ ^{MM}] • 

The domain of L £ + is 7^2 and L+ is a positive self-adjoint invertible operator with a compact 
resolvent. Similarly as for (Li) -1 , we estimate the size of (Li) -1 in C(L 2 ( 



Lemma 2.3 For < e < 1, 

"' r " ,_1 'l£(L2(R)) 



II (£+) 1 ||£(L 2 (R)) ~ e4/3 - 



Proof. See Appendix IA.2[ ■ 
Using Lemma 12.31 we gi ye estimates on various norms of (L^) -1 for sufficiently small e > 0. 
Lemma 2.4 For < e <C 1, 

ll^(^) _1 lk^(R)) < 1, (2-19) 

WULir'WciLHR)) < e 2/3 (2-20) 

\\dx(L £ + )~ 1 \\c(L 2 {R),L^(R)) < e ^ 3 ) (2-21) 

ll(^+)~ 1 |lz;(L 2 (R),L° c (R)) ~ e - (2.22) 

Proof. Let / G L? L 2 and ip = (L^) -1 f . The bound (j2i20l) is obtained by taking an inner 
product of L £ + ip = f with -0 and using Lemma [ 



'\\h(R)+ / 9e|V-| 2 ^<||/|| L2(R) ||^|| i2{R) <e 4 / 3 . 



The bound (|2.22p is a consequence of the bound (|2.20p and Lemma 12.31 applying Sobolev's 
embedding of H 1 (R) into L°°(R) to the function = ^(e 2 / 3 £). To get the bound (pTTUj) . we 
compute 



/ «) 2 dx -2 [ q^'dx + / g 2 ^ 2 dx 
Jr Jr Jr. 

(ip") 2 dx + 2 [ q e (^') 2 dx + [ q 2 ^ 2 dx 



1 > WfWUm = \\L%^\h( 



A p 2 /* 

+ \ i?dx - \ i?dx - -z U 2 {\) + ^ 2 ("1)) , 
£ 2 J\x\<l £ J\x\>l 

where we have used that lirnM^oo q e ipi/j' = 0> which is true because -0 £ D(L £ + ) = Hi- The bound 
flmp holds with the use of the bound (f^22|) and Lemma E3J The bound (j2i2T]) follows from 
Sobolev's embedding of H l (R) into L°°(R) applied to ^'(f) = e 2 / 3 ^'(e 2 / 3 ^) and fr om bounds 
(pT9|) and (|230|) . ■ 
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2.3 The operator (L^) _1 (Li) -1 

Prom the results in the two previous sections, we can deduce easily some estimates on norms of 
(LI) -1 . For instance, 

ll(^+) _1 (^-)~ 1 ||£(L 2 (R)) < ll(^+) _1 ||£(L 2 (R))ll(^-) _1 ||£(L 2 (R)) < ^ ■ 

However, it turns out that these estimates are not sufficient for the proof of the Main Theorem. 
To improve the estimates, we use the fact that if v G B L 2 maximizes [{L £ + )~ l v , v) ~ e 4//3 , 
then (L £ + )~ l v has its L 2 -norm concentrated about the points ±1 (where q £ vanishes), whereas if 
u G B L 2 maximizes ((Ll_) _1 n, uj ~ 1, then (L £ _)~ 1 u has its L 2 -norm concentrated in the interval 
(— 1, 1), away from the points ±1. Figure 1 shows potentials p £ and q e versus x. Figure 2 shows 
schematic shapes of {L £ _)~ l f and (L^) -1 / for a / G L 2 (R). The precise estimates on norms of 
(L^) _1 (L1) _1 are summarized in the following lemma. 




-2-1 1 2-2-1 1 2 

X X 



Figure 1: Profiles of potentials p e (left) and q e (right) versus x. 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 

X X 



Figure 2: Schematic shapes of (L £ _) 1 f and (L £ + ) 1 / for f(x) = exp(— x 2 /4) G L 2 (R). 
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Lemma 2.5 Let a G (0, +00] and 5 > 0. TTien /or < e <C 1, 

ll^)- 1 ^)- 1 !^)) < e 11 / 12 , (2.23) 

IK^r 1 ^) -1 !!^^)) ^ ^ 26/15 ^> (2-24) 

U{\a:\>i}iL%r\Ltr l \\c iL2m < e 7/3 - S , (2.25) 

ll^lxIM-^}^^)" 1 ^-)" 1 !!^ 2 ^)^^)) £ £^(4/3,1/3+30/2)-^ (2 _ 26) 
l|l{|a;|>l-e«} (■£+) 1 (^-) 1 ll£(L 2 (R),L°°(M)) < £ mm(2,l+3a/2)-5 ^ (2.27) 

where if a = +00, we use the convention e a = 0. 

Proof. Let / e%,S= (LIT 1 / and R = {L%)~ l S. We choose 7 G (0, 2/3) (in the sequel, we 
will make different explicit choices of such 7), and we split R into three pieces: R = R\ + R2 +R3, 
where 

R 1 = (L + )- 1 l {W >i } (Li)- 1 /. 

#2 = (^+) _1 l{l-£-r<|2)|<l}(^-) _1 / ; 

R3 = (L e f ) _1 l(_i +e 7 i i_ £ 7)(Li) _1 /- 
Notice that R% and -R3 depend on 7. According to Lemmas 12. 2\ 12.31 and 12.41 

Hflillr3(R) £ ^ 5/3 > ll^ill^W < ^ 7/3 5 H#ilU~(«) < e A/ \ \\R\\\l~w < e 2 - (2.28) 
Thanks to Lemma 12.2^ the Taylor formula provides 

11^11^(1-^,1) £ e 7/2 (l^(l)l + ^II^IIlcC-i,!)) £ e 7/ V /3 + £ ^ (2-29) 
because 7 < 2/3. Thus, using Lemmas 12.31 and 12.41 we obtain 

WR'zhzm < £ 2/3+37/2 , \\R2h m < £ 4/3+37/2 , \\R^ m < ^ 1/3+37/2 , \\m*-<m < ^ 1+37/2 -(2.30) 

The last component i?3 solves the differential equation 

L £ + R 3 = l ( _ 1+£7jl _ £7) S, iel. (2.31) 

We multiply this equality by R3, integrate over R and use the Cauchy-Schwarz inequality. Since 

llS'llraon) i$ 1, we g et 

II4H 2 2 (R) + / Qe\R 3 \ 2 dx < ||i? 3 |k2 ( _ 1+£7)1 _ e7) . (2.32) 
Thus, since ||-R3|li2 ( _ 1+e7il _ e7) < e 2 ~^ f R q £ \R 3 \ 2 dx, 

WR'sWUm + ^ ~ c^ 2 " 7 ) 2 ^ e2 ~ 7 (2-33) 

for some C > 0. We deduce 

WR'sWlHr) < e 1 ^ 2 , \M\v(-i+*r,i-en) Z (2-34) 
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Next, we will establish an estimate on ||-R3||l 2 (i-£t<|2;|<i)- We first estimate the L°°(R) norm of 
R3. Let x De a C°° function on M. with values in [0, 1] such that x( x ) = f° r x < — 1/2 and 
x(x) = 1 for x > 0. We denote the function defined by 



X R 3 (x) := X Rs (l - e 1 + (x - 1 + e^e 1 "^ 2 ) . 



Then, using Sobolev's embedding of H 1 ^— 00, 1 — e 7 ) into L°°(— 00, 1 — e 7 ) (notice that the norm 
of this embedding is the same that the norm of H 1 (R + ) C L°°(M + ), and therefore does not 
depend on e), we obtain 

11-^3 IU°°(0,l-eT) ^ IIX^3||L°°(-oo,l-eT) = 11X^3 ||l°°(-oo,1-£"v) ~ Hff^-cxj^-eT) 

< ^ 1/2+7/4 ||^3||l 2 (-oo,1- £ ,) +£ 1/2 ~ 7/4 ||(X«3) , ||l 2 (-oo,1- £ ,) 

< e 3 / 2 " 37 / 4 . (2.35) 
Similarly, ||-R3||i<x>(_i +£ 7 ) o) < e 3/2 ~ 37 / 4 . Since i?3 solves 

-dlR 3 + q £ R 3 = 0, |x|>l-e 7 , 
where q £ ^ and i?3 G L 2 (R), we infer from the maximum principle that 

P3||l~(r)<£ 3/2 - 37/4 . (2-36) 

On the interval (1 — e 7 , 1), there exists constants C £ A and C £ B such that R% is given by the linear 
combination 

Rs = C%i; £ A + Cm, 
where ip A and ip B are defined in Lemma 12.61 below. 

Lemma 2.6 There exists a constant C > such that for e > sufficiently small, the equation 

- i>"(x) + 2 ^ ~ X =0, -J < x < 1 (2.37) 

/ias two linearly independent solutions ip A and ipg in the form 

^a(x) = a(l - x)Ai (^T^) (1 + Q £ a{x)) , 

= o(l - x)Bi (^^) (1 + , 

where £(x) := ^| Jj^ \/2t(2 — t)dtj , a(x) := (£'(x))~ 1//2 ; Ai ? Bi are i/ie Airy functions, and 
Q A , Q £ B satisfy the bound 

\\Qa\\l°°{-i/2,i) + II<3bIIl o °(-i/2,i) < Ce 2 / 3 . 
Proof. See Appendix IA.3[ ■ 
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According to 10.4.59 and 10.4.63 in [AS], the Airy functions satisfy the following asymptotic 
behaviour at infinity [ASl Section 10.4]: 

Ai(2) ~ ^tW" § " 3/2 and Bi(z) ~ ^h^ e¥m as ^+°°- ( 2 - 38 ) 

At the point x = 1, we deduce from (|2.36p that 

|Cio(0)Ai(0)(l + Q%(1)) + C£a(0)Bi(0)(l + Q%(1))\ < e 3 / 2 ' 3 ^. 

Thus, 

Km I < e 3/2 -'^ /4 + \C%\. (2.39) 
At the point x = 1 — e 7 , provided that 7 < 2/3, we similarly have 



C%a(e^)M ( ^1) (1 + Q\{1 - e 7 )) + C%a(e^)Bi (*&± ) (1 + 1; ( 1 : ) 1 



< r 3/2-3 7 /4 



Since 

~ 2 2/3 x as s (2.40) 
and thanks to (|2.38p and (|2.39p . we obtain 



£ 3/2-3 7 /4 



\Cb\Z-, r and \C%\ < e 3/2 ~^ /4 , (2.41) 

where Bi ^^7^ — > 00 as e — > 0. Since 7 < 2/3, one can choose € (7, 1 — 7/2). Using again 
the maximum principle, we get 

\R 3 (x)\ ^ \R 3 (1 - £i + xyl-e^ + e 13 . 

Moreover, thanks to (|2.4ip . we have 

\RM - £ 1 + e /3)| < e 3/ 2 -37/4 Ai f gfc 7 -^ 3/2-3 7 /4 Bl ( Q 



Using (I2.40p again, we deduce from (I2.38|) that there exist a constant cq > such that 
£ 3 / 2 ^/ 4 Ai ( a£ l~f ] ) < exp (-coe^ 2 - 1 ) , 



Bi 



£ 3/2-3 7 /4_V^^i < (_ /3+ 7 /2-l 
Bi ^ 



where we have used 
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which holds because (3 £ (7, 1 — 7/2). Therefore, we find 

II^IL-a-^.+oo) £ 1^3(1 - h 7 + eP)\ < exp (-coe^^ 2 - 1 ) , (2.42) 

which shows that ^3(1) and C\ are actually exponentially decaying as e — > 0. Then, we infer 
from (12361) and (I2~l2"1) 

11-^3 ||i2(l- £ 7,l) < ll^3|lz / 2( 1 _ e7>1 _ £7+£/ 3) + ||-R3|| L 2( 1 _ e7+£/ 3 i1 ) 

< e/ 3/2 £ 3/2-3 7 /4 + £7 /2 gxp ^^^+^2-1^ < £ 3/2 +/ 3/2-3 7 /4_ ^ 

The L 2 norm of R3 on the interval (— 1, — 1 + e 7 ) is estimated in the same way. Next, we estimate 
the I? norm of -R3 on the interval (1, 00). We multiply (|2.31|) by -R3 and integrate over (1, +00). 
Since p £ ^ 1 for x ^ 2 and e ^ 1, we obtain 

/+00 r+00 
(R' 3 ) 2 dx + J^ p £ Rldx = -R 3 (l)R' 3 (l) 

< expf-co^/ 2 - 1 )^ 3 , (2.44) 



where -^3(1) has been estimated with (|2.42p and the bound for R 3 (l) comes from Lemmas [27 
and[2J~J The I? norm of R 3 on (1,2) is estimated thanks to (12.42H . Together with (j2.44[) . we 
deduce that 

ll^3lk 2 (i, + oo) < exp (-c^+^Z 2 - 1 ) , 
where c = co/2. The 1? norm of R 3 on (—00, —1) is estimated similarly, thus 

\\R3\\m\ x \>i) < exp (-ce^/ 2 - 1 ) . (2.45) 

Since R 3 solves 

-R'i + q £ R 3 = 

on (1 — e 7 , +00) and -R3 £ L 2 (R), we deduce from the maximum principle that if -R3 does not 
identically vanish on (1 — e 7 , +00), then R 3 has a constant sign on that interval. For instance, 
i?3 > (the argument is similar in the other case). Then, R 3 (x) ^ for every x ^ 1 — e 7 . 
Therefore R 3 is a negative increasing function on (1 — e 7 , +00). Let us assume by contradiction 
that \R' 3 (1 - e~* + > exp(-c e /3+7/2 ~ 1 )/e 2 . Then, for x ^ 0, it follows from the Taylor 
formula and (j2.42j) that for e sufficiently small, 

i? 3 (l-e 7 + e /3 + e) = «s(l - + e") + eR' 3 (l - e 7 + e p ) + / / R' 3 \l - e 7 + e p + t)dtds 



R 3 (l - e 7 + e 13 ) + eR' 3 (l - e 7 + e p ) + f f R' 3 \l -e^ + e^ + t) 

Jo Jo 

,(-0^+7/2-1) f C 



^ exp(-c ^ +7 /- 1 ) ( C - ^ + C £ -^) < 0, 
for some C > 0, which is a contradiction with the positiveness of R 3 . As a result, 



|-^3ll-L°°(l-eT+e' 3 ,+oo) 



R' 3 {l-e< + eP) <exp(-ce /3+7 / 2 - 1 ). (2.46) 
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At this stage, we have established all the estimates required to prove the lemma. First, (|2.28|) . 
(123(1 and (12341 yield 

\\ R lli 2 (K) ^ II-R'iIIl 2 (R) + \\ R 2 IIl 2 (R) + \\ R 3 \\l 2 (R) 

< £ 5/3 + £ 2/3+3 7 /2 +£ l- 7 /2_ (2 _ 47) 

The choice 7=1/6 provides (12331) . From (1235]) . (12301) . (12341) . (12331) and (123SD . we obtain 

II-R||l 2 (R) < H- R i||l 2 (R) + H- R 2||l 2 (R) + l|- R 3||L 2 (-l+e7,l- e 7) + ||#3||l 2 (1-£7<|;c|<1) + 11-^3 IIl 2 (|x|>1) 
< £ 7/3 + £ 4/3+3 7 /2 + £ 2- 7 + £ 3/2-3 7 /4 +/ 3/2 + gxp ^_ c£ /3+ 7 /2-l^ _ (2 _ 4g) 

The choice 7 = 4/15,/? = 13/15 — 25, for sufficiently small positive number 5, provides the bound 
(I2.24|) . Similarly, we have 

l|-R||i 2 (|a:|>l) ^ ll^l|li 2 (|x|>l) + 11-^2 |U 2 (|a:|>l) + 11-^3 lli 2 (|x|>l) 
< £ 7/3 + £ 4/3+3 7 /2 +exp ^_ c£ /3+ T /2-l^ 

The choice 7 = 2(1 — 5)/3, (3 = 2/3, for any small positive number 5, provides the bound (|2.25|) . 
If a > 0, 7 < min(a,2/3) and if e is sufficiently small, we also obtain from (|2.28p . (12.30P and 

HMD, 

ll-R \\L ao (l-e a ,+oo) ^ \\$\\L°°(l-ei+eP,+oo) ^ I|-RiI|l°°(R) + II^IU 00 ^) + 11-^3 Wh^il-e^+el 3 , +00) 

< £ 4/3 +£ l/3+3 7 /2 +exp ^_ c£ /3 +7 /2-l^ {25Q) 

A similar argument on (— 00,— 1 + e a ) gives (|2.26p . for the choice 7 = min(a,2/3) — 25/3, 
f3 = (1 +7)/4. If 7 < min(a,2/3), thanks to ([535]) . ([230]) . ([2351) and its twin estimate on 
(—00, —1 + e a ), we get similarly, for e sufficiently small, 

ll^llL°°(|x|>l-e a ) ^ H-R|lL°°(|a;|>l-e^+£' 3 ) ~ 1 1 -^1 1 1 L°° (R) + 11-^2 ||l°°(R) + II ^3 1 1 (|a;| >l- £ -K+e^) 

< £ 2 + £ l+3 7 /2 + exp ^_ Co£ /3 +7 /2-l^ (2 51) 

The bound ([237]) follows from ([531]) . again with the choice 7 = min(a, 2/3)— 2<5/3, /3 = (l+7)/4. 



3 Proof of the Main Theorem 

3.1 The operator A e for e > 

We consider here the operator 

A e := e~\-dl+p £ (x))-\-dl + q^x))- 1 = e~ 2 (Ll)" 1 (L^y 1 . (3.1) 

As we have seen before, if e > 0, both operators LI and L £ + on L 2 (M) are invertible with compact 
resolvent. As a result, A £ is a compact operator on L 2 (M) for any fixed e > 0. Thus, its spectrum 
consists of a sequence of eigenvalues which converges to zero. Moreover, these eigenvalues are 
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all strictly positive. Indeed, if fi is an eigenvalue of A £ and if is an associated eigenvector, 
C := (L £ + )~ l / 2 ip satisfies 

(L%)- i / 2 (L £ _r i (L%r 1/2 c = K. 

Therefore, /i is an eigenvalue of the self adjoint positive operator (L £ ^) _1 / 2 (Li) _1 (L £ h ) -1 / 2 , which 
implies /i > 0. We order eigenvalues of A e as 

()<•■■< Hn,e ^ ■ • ^ (J-2,6 ^ Ml.eps < °0. 

3.2 The operator A 

As e — > 0, we can formally expect that A £ converges in some sense to the operator 

A = (-dl+p r 1 ' 



2(1 -x 2 Y 
where 

if|x|<l, 
+oo if \x\ > 1. 

Let us describe more precisely the action of the operator Aq on L 2 (R). The following lemma is 
helpful for that purpose. 



Po(x) = 



Lemma 3.1 If u € L 2 (R), then (j^p) G (# 2 n #<})' (-1, 1), w/iere (# 2 n H%) (-1, 1) 



endowed with the H 2 norm. Moreover, the map u ^ (jr^?) ^ s continuous from L 2 (1R) into 

Proof. By Sobolev's embedding theorem, H 2 (—l, 1) is continuously embedded into C 1 ([— 1, 1]). 
Therefore, if g € (H 2 n #o)(-l, 1), then 

\g(x)\ = \g(x)-g(±l)\^\\g'\\ L -(l-\x\), 

with +1 for x > and —1 for x < 0. It follows that for every x € (—1, 1), 



5^) 



1-x 2 



^ - — — ^ IMIfra- 



1 + x 



As a result, using the Cauchy-Schwarz inequality, we obtain 

tt(x) 



1 — x : 



-g{x)dx 



< 



\ U \\L 2 (R) 1 1 PI Iff 2 (—1,1) > 



which completes the proof. 



Let us denote the Dirichlet realization of the Laplacian A = d 2 on the interval (—1,1) by Ap. 
It is well known that (— A^) -1 maps continuously L 2 (— 1, 1) into (if 2 n -ffg) ( — lj !)• By duality, 
it also continuously maps (H 2 H H^)' (-1,1) into L 2 (-l, 1). For u £ L 2 (R), A u £ L 2 (R) is 
defined by 



(A u)\ { \ x \ >1} = 0, 

;(^),(-i,i)=(-a d) -((^ t ) i( _j 



(3.2) 
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Thanks to Lemma O and the continuity of (-Ad)" 1 : (H 2 n H$)' (-1, 1) i-> L 2 (-l, 1), ^ is a 
bounded operator on L 2 (M). Moreover, we have the following lemma. 

Lemma 3.2 For any u G L 2 (R) and any s£ [-1,1], 

• 4o " (<,) = /, {Lw^^-l wr*n iv+ — I{u) < (3 - 3) 



where 



In particular, Aqu is continuous on R. 
Proof. For any u G L 2 (R) and any y G (—1,1], we have 
^ u(x) 



{l + x 



-dx 



1 1/2 , „i 1 v 1/2 II || 

1 1 ' Wi) H U f? 4 ) « W ■ (3 - 5) 



u 



which implies that the map u \— > J^ 1 y-^dx is continuous from L 2 (R) into L x (— 1, 1). Similarly, 

one can see that the map u ^> f\ ^z^dx has the same property. As a result, n i— > I(u) is 
a continuous linear form on L 2 (R), and the map which assigns to u the right hand side in 
(13. 3p is continuous from L 2 (R) into L°°(— 1,1) C L 2 (— 1,1). As we have seen before, so is 
u i— ► (^o w )|(-i,i)- Actually, both sides in f|3.3[) only depend on the restriction of u to (—1, 1), so 
that they can be considered as continuous from L 2 (— 1, 1) into itself. Therefore, using the principle 
of extension for uniformly continuous functions, it suffices to check (|3.3p for u in a dense subset 

ofL 2 (-l,l). This can be done for u G C£°(— 1, 1). Indeed, in this case (jr^z) ^ G L 2 (-l, 1), 
therefore (A) u )|(-i,i) G (-f^ 2 H -ffg) ( — 1, 1)- In particular, lim ^(Aoit)(s) = 0. On the other 
side, we can easily check that the right hand side in (|3.3|) also vanishes at s = ±1 and its 
second derivative is — 2 (i^i) > which completes the proof of (|3.3|) . It remains to prove that 
lim (Aqu)(s) = is true for any u G L 2 (R). This follows from the fact that the maps 

s^±l=F0 

y ^ ly T+^dx and y h-> j^dx are in L l (-l, 1). ■ 



Lemma 3.3 Aq is a compact operator on L' 

Proof. By Lemma 13.21 Aq is continuous. Thus, according to a standard criterion of relative 
compactness for a subset of L 2 (R) (see, for instance, Corollary IV. 26 in [B]), it is sufficient to 
check the following two conditions 

(i) for every n > 0, there exists a compact subset u> C R such that for every u G B L 2, 

\\ a ou\\l 2 (r\u) < V 

(ii) for every rj > and for every compact subset wCl, there exists 5 > such that for every 
u G B L 2 and for every h with \h\ < 5, 

\\A u(- + h)- A u\\ L 2^ < n. 
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In our case, condition (i) is trivially satisfied: we choose u = [—1,1] and then ||A)^||l 2 (]r\(j) 
for every u G B L 2. To check condition (ii), we note that if— 1 sj s,s + h ^ 1, then 



\A u(s + h) - A u(s)\ 



/ / 77 - / ' dx dy + -I(u) 

J s Vi i 4(1 -x) J v 4(1 + x) J 2 y ' 



s+h 



< V\h\.C\h\ 



u \\L 2 (R) ( 1 1 



+ 75 dy 



\h\ 



2 c|I u IIl2( 



4 2' 

for some constant C > 0. A similar estimate holds if either +1 or —1 lies between s and s + h 
(which can only happen if \s\ < 1 + \h\), whereas if both s and s + h are outside of (—1, 1), then 
A u(s + h) — Aqu(s) = 0. Therefore, 



\\A u(- + h)- A u\\ L2m < (2(1 + \h\)fl 2 I vE. + j , 

and condition (ii) follows. ■ 

Since Aq is compact, its spectrum is purely discrete. Clearly, is an eigenvalue of Aq and 
the associated infinite-dimensional eigenspace is made of the set of functions in L 2 (R) supported 
in the exterior domain {x E R : \x\ ^ 1}. If fi ^ is an eigenvalue of Aq and w an associated 
eigenvector, it follows from the definition of Aq that w = on {x G R : |x| ^ 1}, whereas on 
{x E R : |x| < 1}, w solves 

- 2(1 - x 2 )w"(x) ='jw(x), -1<x<1, (3.6) 

where 7 = 1/fJ,. Moreover, thanks to Lemma 13.21 w = jAqw is continuous so that w(—l) = 
w(l) = 0. We shall now prove that the only solutions of (|3.6p vanishing at the endpoints ±1 are 

— 1/2 

the Gegenbauer polynomials C n+1 (x) for 7 n = 2n(n + 1), where n ^ 1 is integer. Thus, the 
spectrum of operator Aq is given by 

a(A Q ) = {fin := 7 1 —r , n^llu{0}. 
[ 2n[n + 1) J 

1/2 

Lemma 3.4 Equation S3. 6\) admits a family of solutions ("f,w) = (j n ,C n+1 ), forn ^ —1, where 

7„ = 2n(n+l) and is a Gegenbauer polynomial with degree m. If (7,10) G" {(7 n , aC n M 2 ) | n ^ 
— l,a G R} is a solution of S3. 6\) . then it satisfies 

lim (\w(x)\ + \w(-x)\) ^ 0, lim (\w'(x)\ + \w'(-x)\) = 00. (3.7) 

:c^l— z^l— 

The only solutions (7, if) 0/ (f ff. 6p suc/i i/iai w(l) = w(— 1) = are (7 n ,aC n+1 ), /or n ^ 1 and 
a G R. 

— 1/2 

Proof. Explicit computations show that Gegenbauer polynomials C n+1 (x) from Section 8.93 
in |GR] are solutions of (|3.6p for 7 n , for any n ^ — 1. In particular, for n ^ 1, by equation 8.935 
in [GRj . we have 
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which proves that C n V 2 (l) = C n V 2 {— 1) = for n ) 1, whereas C 1 ^ 2 (x) = 1 and C l 1 ^ 2 {x) = 

_ ]/2 

—x. We next prove that if (7, w) solves (|3.6p and w is not proportional to C n+1 with n ^ — 1, then 
w satisfies (13. . We introduce the new variable z = x 2 for < x < 1, and the function -u(z) := 
tu(x). It is equivalent for w(x) to solve ()3.6[) on (0,1) or for -u(z) to solve the hypergeometric 
equation: 

z(l - z)u"(z) + -(1 - z)u'(z) + 1u{z) = 0, < z < 1. (3.8) 
2 8 

This equation admits a general solution given by 9.152 in |GR| 

u(z) = Cl F(a, b, c; z) + c 2 z 1/2 F (a + ^ b + ^ ^; z^ , (3.9) 

where 

a + o= — , ab = , c=- 

2 8 2 

and F(a,b,c; z) is a hypergeometric function. Clearly, the function x 1— » u(x 2 ) = -u;(x) defined 
by (|3.9p is analytic for < x < 1 and can be extended into an function w which is analytic for 
— 1 < x < 1, given by 

/ l 1 3 r 

w(x) := aF(a,b,c;x 2 ) + c 2 xF [ a + -,b + -,-;x' 



2' 2'2 

where the first term is even in x and the second term is odd in x. Since w solves (|3.8p . the 
uniqueness in the Cauchy-Lipshitz Theorem ensures that w = w. In order to prove the Lemma, 
it is sufficient to consider one component of the solution at one boundary point, e.g. F(a, b, c; x 2 ) 
at x = 1 (z = 1). Since Re(c — a — b) = 1 > 0, the function F(a, b, c; z), which is analytic on 
{z : \z\ < 1}, is also bounded as z — > 1 (see 15.1.1 in |AS| ). Using 15.1.20 in [AS], that is 

, , T(c)T(c-a-b) 

F(a,b,c;l )= )> \ -(, 

T(c — a)T(c — b) 

we find that 

tt 1 / 2 sin(7ra)r(-a) _ cos(vra)r(l/2 + a) 

(a ' ' q } ~ r(l + a)r(l/2-a) ~ ~vr 1 /2r(l/2-a) ~ vr 1 /2 r (l + a) ' 

Parameters a and 7 are related by 7 = 4o(l + 2a). If 7 = 72m-l = 4m(2m — 1) for m ^ 1, then 
either o = — m or a = —1/2 + m, both give F(a, b, c; 1) = 0, corresponding to even polynomial 

— 1/2 

solutions C 2m ■ For all other values of 7 and a, F(a, 6, c; 1) is bounded but non-zero. On the 
other hand, using 15.2.1 in [ASj . that is 

-^-F(a, b, c; z) = —F(a + 1, b + 1, c + 1; z), 

az c 

since Re(c + 1 — a — 1 — 6 — 1) = 0, we obtain that 4;F(a, b, c; z) = 2x-^F(a, b, c; z) diverges 
as z — > 1 (see 15.1.1 in [AS] ), unless the series for F(a,b,c,z) is truncated into a polynomial 
function, which happens precisely when a or b is a negative integer, which implies that 7 equals 
one of the 72m-i's for some m ^ 0. Therefore, lim^^i |-u/(x)| = 00 if w(x) is an even solution of 
(13. 6p and 7 7^ 72m- 1 f° r m ^ 0- Similarly, the statement is proved for an odd solution of (|3.6p . 
given by si 7 (a + 1/2, 6 + 1/2, 3/2; x 2 ) for 7 7^ 72™ with m ^ 0, where 7 = 72m = 4m(2m + 1) 

— 1/2 

correspond to odd polynomial solutions C 2m _ 1 . ■ 
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3.3 Convergence in norm of A £ to A as e — > 

Our goal in this section is to prove the following result. 
Theorem 3.5 It is true that 

A £ — ► A in C(L 2 (R)) as e -> 0. 
Once this result has been proved, we immediately have the corollary. 
Corollary 3.6 For every integer n ^ 1, 

Vn,e — ► Mn as e ^ 0. 

Moreover, if w n is an eigenvector of Aq associated to the eigenvalue /%, there exists a set 
( w n,e)e>o C L 2 (M) of eigenvectors of A £ associated to the eigenvalues fj, n£ for e > 0, such that 

w-n,e — ► w n in L 2 (M.) as e —> 0. 

Proof. Since convergence in norm in C(L 2 ) implies generalized convergence, it follows from 
Theorem 3.16 on p. 212 in [K] that for every integer N ^ 1 and for < e <C 1, 

Mat + VN+i . \ _ , A 
, +oo Pi a\A £ 



- N. 

Moreover, fi n£ — > fi n as e — > 0, for any 1 ^ n ^ N, which proves the convergence of the 
eigenvalues. For the eigenvectors, let us fix n ^ 1, and let fi„ C C be a neighborhood of \x n 
such that Q n does not contain nor any other eigenvalue of Aq. From the convergence of the 
eigenvalues, it follows that for e sufficiently small, A £ has a unique eigenvalue in fi n , which is /% )£ . 
For any integer m ^ 1, we denote by E m (resp. E^) the eigenspace of Aq (resp A £ ) associated 
to the eigenvalue [i m (resp /%i )£ ). We also define 

F n := ( © E m ) 0Ker^o and F rit£ : = © E £ m , 

as well as P n 6 £(L 2 (M)) (resp. P n ,e) the projector on E n (resp S n ,e) along F n (resp. F nj£ ). 
Then, Theorem 3.16 in [K] also ensures that P n>£ — ► P n in C(L 2 ) as e — > 0. Thus, w nfi := P n ^ £ w n 
is an eigenvector of ^4 £ for the eigenvalue (J, n , £ , and we have 



(Pn,e - Pn)Wn\\L 2 (R) ^ ||-Fn,e _ Pn\\c(L 2 (S.)) \\ w n\\l?(M) °, 



which completes the proof. 



Remark 3.7 A straightforward consequence of Theorem \3.5\ is that A* — ► Aq in £(L 2 (1R)) as 
e — > 0. Thus, an analogous result to Corollary \3. 6\ holds for the eigenvalues and eigenvectors of 
At and An 



1 o- 
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The convergence statement of the Main Theorem directly follows from Corollary 13.61 since 
the spectrum of system (jl.5p is made is made of the eigenvalues A = ±ie/y/Jl, where \i describes 
the spectrum cr(A e ) of A £ . Indeed, if (A, u,w) G C x L 2 (R) x L 2 (R) solves (|1.5j) , a straightforward 
computation shows that 



A e w = -j2 w ' 

thus A = for some \i G a(A e ). Conversely, if A £ w = [aw, with w G L 2 (M), then (ie/^/Jl, u, w) G 
C x L 2 (R) x L 2 (M) solves system (TT3]) with 

Let us now turn to the proof of Theorem 13.51 In order to compare Aqu and A e u for e > 
and it G L 2 (R), we would like first to express Aqu as Aqu = A £ (A £ )~ 1 Aqu. This can be done 
with the help of the following lemma. 

Lemma 3.8 Let H be a Hilbert space and L be a self-adjoint operator on H with domain D{L) 

endowed with the graph-norm \\ ■ \\d(L) = (II " \\% + 11-^ " IIh) 1//2 - Assume that L is continuously 
invertible and X is a Banach space continuously embedded in H. L induces an operator Lx on 
X, defined by 

D(Lx) = {x G X, Lxx G X}, Lxx = Lx for any x G D(Lx)- 

1 /2 

D[Lx) is endowed with the graph-norm || • Htv^) = (|| • \\\ + ll-^x " llx) ■ Assume further that 
D(Lx) is dense in H and that D{L) is continuously embedded in X . Then L is extended to X' 
as a bicontinuous map Lx> ■ X' ^> D(Lx)' defined by 

(Lx<f,<p) D ( Lx y iD ( Lx ) ■= (f,L X (p) X ',x f or an V f ^ X' and if e D(L X ). 
Proof. See Appendix IA.41 ■ 

To prove that Aqu = A £ (A £ )~ 1 Aqu for any e > and u G L 2 (M), we apply Lemma 13.81 twice. 
For the first application, H = X = L 2 (1R) and L = L £ _, such that L £ _ is extended as a bicontinuous 
map (also denoted L £ _ for convenience) from L 2 (M) into D(L e _)' . Thus, Aqu = (L e _)~ l L £ _Aqu. 
For the second application, H = L 2 (R), X = D(L £ _) and L = L £ + such that L £ + is extended as a 
bicontinuous map (that we will also denote L £ + ) from D{L £ _)' into 

D D{ L t )(L £ + ) := {v G D(Ll), L%v G D(L £ _)}. 

Note here that D(L £ + ) is continuously embedded in X = D(L £ _), since L^ — Lt = -^-pr^l(_i i) G 
C(L 2 (R)) (actually, D(L £ + ) = D{L £ _) and the norms || • \\d(l e ) an d || ■ ||z>(l s ) are equivalent). As 
a result, 

A Q u = {L £ _)-\L £ + )- 1 L £ + L £ _Aqu = A £ e 2 L £ + L £ _A u = A £ (A e )- x A Q u, 
where (A £ )~ l maps D D ( L ? )(L £ + ) into L 2 ' 
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The identity (|3.3|) provides an explicit expression of Aqu for any u G L 2 (R). Let us next use 
this identity to express L £ _Aqu G D(L £ _)' . If <£> G D(L £ _) and it G L 2 (M), then direct computations 
involving integration by parts give 



(L £ _A u^) DiLtyMLt) 

= (A Q u, Ltip) L 2 L 2 = - 
- 1 u(x) 



(A u)(s)ip"(s)ds 



4(1 + z 



1 w(x) 
4(1 + x 



■dx 



u{x) 
-1 4(1 -x) 
u(x) 

! 4(1 - X 



I + 



s- 1 



I (it) ) (p"(s)ds 



dx I tp'(s)ds 



I(u) 



Ml)-^(-l)). (3.10) 



Performing another integration by parts, the first term in the right hand side of (|3.10|) can be 
expressed as 



u x 



■dx 



u(x) 



lim 

5^0 



lim 

<5-+0 



1+5 
1 



1 u{x) 



4(1 + x 

1-8 

4(1 + x 
u{x) 
l+s 4(1 + x) 
1 u(x) 
1 4(1 + x) 



i4(l-x) 



dx I ip'(s)ds 



■dx 



5 u(x) 
1 4(1 - x) 



dx I tp'(s)ds 



(ip(x) - <p(-l + S))dx + 



1-6 



u(x) 
1 4(1 - x] 



((f(x) — ip(l — 6)) dx 



(ip(x) - <£>(— 1)) dx + 



1 u(x) 
1 4(1 - x] 



(<p(x) -<p(l))dx. 



(3.11) 



The first limit in the right hand side of (|3.1ip is evaluated as follows. (The second limit is 
evaluated similarly.) We write 



u(x) 



_ 1+s 4(1 + x) 

1 u(x) 



(</?(x) — tp(— 1 + 5)) dx 



1 u(x) 

.1 4(1 + x; 



(</>(x) -ip(-l))dx 



I 7#tM-i)-^-i + *))^- / 

7-1+5 4(1 + x) J_ x 



l+<5 



it(x) 
4(1 + x) 



(<p(x) - ¥?(-!)) dx 



.(3.12) 



The two terms in the right hand side of (|3.12p converge to as 5 goes to thanks to Lebesgue's 
dominated convergence theorem. For the first term, the integrand is dominated by 



u(x) 



4(l + a?; 



M-l)-^(-l + 5))l ( _ 1+5jl) 



5u{x)\\<p'\\l°° 
4(1 + x) 

■u{x)\\W\\ L ™ 



GL^-l,!). 



The integrand of the second term is dominated by the same integrable majorant. Then, from 
(PLIOJ) and (gEHD we deduce that 



[L E _A U, <P) D ( L e_y tD ( L eJ 

^ u(x) ip(x) — <p(—l) 



(3.13) 



1 + x 



dx + 



1 4 1 — x 2 
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Thus, if u G L 2 (R) and (p G D D(Li) (L^), then 



),D(L S _) 



u(x)ip(x)dx 



u(x) 



l4(l + x) 



+ 



(^"(1) -</(-!)) 



1 

-i 4(1 -x) 
u(x)p(x)dx. 



<p"{x) - <//'(!)) dx 



\x\>l 



Finally, if we introduce the adjoint operator of A e , 

A* £ := ^(LirHLtr 1 G C(L\R),D D(Lt) (L%)), 

we get for any u, <p G L 2 (M) 

(A u - A E u,(p) L 2 iL 2 = (A E (e 2 L £ + L £ _A u- u),p) L2 ^ 2 
= (e 2 L £ + L £ _A u - u,A^) DD(Ltmy>DD(Le _ m) 

_ £ 2 f 1 u(x) (A* £ cp)"(x)-(A^)"(-l) dx g2 f 1 n(x) (A^r(x)-(^^r(l) ^ 



+ 



i 4 
e 2 I{u) 



1 + x 

{(A^)"{1) - (A^)"(-l)) 



\x\>\ 



-i 4 



1 - X 



(3.14) 



In order to prove the convergence of A £ to Aq in £(L 2 (]R)), it is sufficient to prove that the right 
hand side in (|3. 14[) converges to as e — > uniformly for u,<p G B L 2 . Up to terms which may be 
estimated similarly, it hence suffices to prove that the three quantities 

Q\(u,v) := \e 2 I{u){A* £ p)"{l)\, 



Q £ 2 (u,ip) :-- 



\x\>l 



u{x){A*p){x)dx 



£2 / 1 (A^M-iA-vrm^ 

1 — x 



defined for u^e L 2 (R), converge to as e — > 0, uniformly for u, (p G B^2. In other words, we 
should choose u and p in B L 2 and prove that 



Q\{u, <p) + Q E 2 (u, ip) + Q e 3 (u, if) < C(e), 
where C(e) does not depend on u or <p and C{e) — > as e — > 0. 



(3.15) 
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Estimate on Q\ . We have already seen in the proof of Lemma 
other side, 



that \I{u)\ < 1. On the 



e 2 ^ = g B (^)- 1 (Li)- 1 -(iii)- 1 . 
Since <fe(l) = 0, it follows from Lemma 12.21 that 

|(^flg^)(l)| = |((tf-)~V) (l)|<e 2/3 . 

Therefore 



(3.16) 



Estimate on Q\. It follows from Lemma 12.51 and from the Cauchy-Schwarz inequality that 

Q§(«,¥>) < ^ /3 ~ S , (3-17) 

for any 5 > 0. 

Estimate on Q%. Thanks to the Cauchy-Schwarz inequality, it suffices to prove that 



(e 2 d 2 A*Mx) - (e 2 d 2 A* £ Ml) 



1 — x 



as e —* 0, 



L2(-l,l) 

uniformly for ip £ -B^2. Using a commutator, we first decompose the operator e 2 l(_i^d 2 A* as 

e'l^dlAt = -l { _ X!l) Ll{L%)-\LL)- 1 

-l^L'+y 1 + l^d 2 [{L%)-\ (Li)- 1 ] . (3.18) 

We introduce the functions r := (L e + )~ l ip, s := (Ll) -1 ^, R := (L^) _1 s, S := (L £ _) r and 
w := d 2 {R- S). Then, 



(e 2 c^M*) - (e 2 <9^Ml) 



1 — x 



L2(-l,l) 



r(x) — r(l) 



1 — x 



+ 



L2(-l,l) 



cj(x) — 



1 — X 



L»(-l,l) 



According to Lemma |2.4| ||r'||^oom) < e 1 ' 3 and the first term is hence estimated by 



r(x) — r(l) 



1 — x 



< pV3 



(3.19) 



Let us now estimate the second term in the inequality above. If we make the difference of the 
two fourth-order differential equations satisfied by R and S on (—1, 1), we find that to solves the 
differential equation 



2(1 -x 2 ) 



-co 



—^R H — ttR'i 



■1 < x < 1. 



(3.20) 



Let a £ (0, 2) (different explicit choices of a will be made later), (3 = 23/30 — 5 and 7 = 7/15 + 6, 
where < 5 < 1/45. Thanks to the triangle inequality, 



Uj(x) — U>(1) 



1 — X 



< 



Il 2 (-i,o)+ £ 7 II w IIl2(o,i- E 7) +e 7 |w(l)| + 



w(:e) — U>(1) 



1 — X 



(3.21) 

L 2 (1-£T,1) 
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Next, for x G (— 1, 1), we have 



u(x) =d 2 x {R- S){x) = r{x) - s(x) + 2(1 fi- R(x) (3.22) 



and 



J(x) = r'ix) - s'(x) + 2(1 fi- R'(x) - %R{x). (3.23) 
Thanks to Lemmas 12. 2\ 12. 4| and 12.51 we obtain 

|u;(±l)| = |r(±l)-s(±l)|<e 2 / 3 (3.24) 

and 

|o/(±l)| 



r'(±l) - s'(±l) T J^(±l) 



<l + ^< £ - (3.25) 



If we multiply (|3.20p by u, integrate over (—1, 1) and use the Cauchy-Schwarz inequality, we get 

1 I' 1 

IkllV!,!) + ^ J (l-X^dx 

^ ^2 ll B IU 2 (-i,i)IMU 2 (-i,i) + ^II^IIl 2 (-i,i)IMU 2 (-i,i) 

(1)1 + |w(-l)||u; (-1)| + j 2 . (3.26) 

Decomposing (—1, 1) into (—1 + e a , 1 — e a ), (—1, —1 + e a ) and (1 — e a , 1) and using the Taylor 
formula and the Cauchy-Schwarz inequality on the last two intervals, we get thanks to (|3.24p 

IMIl 2 (-i,i) < |M| i2( _ 1+£ « + e a ' 2 (|w(l)| + |w(-l)| + ^llo/llia^!)) 

< ||o;|| x2( _ 1+£ « + £ Q / 2+2 / 3 + e a ||a/|| La( _i i i ) . (3.27) 
Prom (ET26]) . (I3T241 . (I3T25]) . flg]gTP and Lemma [23] we deduce, for sufficiently small 6 > 0, 

II, /II 2 1| i|2 

ll w IIl 2 (-1,1) + e IMIl^-l+e^l-e") 

< ^26/15—5—2 / II, ,|i , ^a/2+2/3 , _ai|, ,/|| ill, /II \ , ^2/3-5 , c .2/3-<5 
~ £ (JMIl 2 (-1+£<*,1-£ q ) + e + £ IMIl 2 (-1,1) + IMIl 2 (-1,1)J + e + e 

< i _a/2+2/5-(J , _-4/15-5||, ,11 , _-4/15-5 II. /II /a oq\ 

Therefore there exists a positive constant C such that 



l w IIl 2 (-i,i) - Ce I +£ I IMlL 2 (-l+£«,l-£ a ) ~ Ge 



< g 2 / 3 " 5 _|_ £ a/2+2/5-<5 _|_ g -8/15-2<5 _|_ £ 22/15-a-2<5 (g 29) 

We deduce that for any a £ (0, 2), 

< £ 26/15- a - S + £ 4/3-«/2-8/2 + e 6/5~ a /4-5/2 + £ ll/l 5 -a/2-8 

< e "/i5-«/2-i ( 3.3 0) 
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and 



W\\ti><.-l,l) £ £- 4 / 15 - 5 + e V3-« + £ l/5+a/4-5/2 + £ -4/15-« + £ ll/5-«/2-« 
< -4/15-5 

Using d32Z])> dSMD , and (JMH), we obtain 

II, ,|| < _ll/15-a/2-<5 t _Of/2+2/3 i _-4/15+a-<5 

IMIl 2 (-1,1) ^ e + e 

For a = 2/3, we get 

II, ,ll < =3/5-1$ 

IMIl 2 (-i,i) ^ e 

Coming back to (|3~2T]) . thanks to (g^H) , (|33UD with a = 7, and (ET32"|) . we obtain 



(3.31) 



oj(x) — cj(1) 



1 — X 



< e 2/5-5 + £ ll/15-3 7 /2-<5 + £ 2/3- 7 + 

L2(-l,l) 

If 7 = 7/15 + <5 and (3 = 23/30 - 5, we have 

1 _ e 7/15 + £ 23/30-5 < ! _ £ 7 

for sufficiently small e > and therefore 
u;(x) — to(l) 



1 — x 



L 2 (1- £ V) 



(3.32) 
(3.33) 



1 — x 



u(x) — cj(1) 



1 — X 



L 2 (l-eV) 

From (|3.22p we infer, for x E (—1, 1), 

w(x) - _ r(x) - r(l) s(x) - s(l) i 2(1 + a;) 
1 — x 1 — X 1 — X 

Like in (|3.19p , it follows from Lemmas 12.21 and 12.41 that 

r(x) — r(l) 



L 2( 1 _ £ 7/15 +£ 23/30-«5 il ) 
, 12(x). 



+ 



1 — X 

s(x) — s(l) 



< £ 17/30 



L 2(- 1 _ £ 7/15 +£ 23/30-5 jl ) 



1 - X 



< P 7/30 



(3.34) 

(3.35) 
(3.36) 



L2( 1 _ £ r/15 +£ 23/30-« >1 ) 

Splitting i? as Ri + R 2 + R3 as in the proof of Lemma [231 and using ([2T28]) . (pT30|) and 
we deduce that 



1 1211 



| i2(1 _ £ 7/i 6+£23 /30- 5il) < £ 7 / 3 + £ 61 /30 + 35/2 + exp (_ c£ 23/30- 5+ 7/30-l) < £ 61/30 j (3.37) 

for some c> 0, since 7/15 < 2/3 and 7/15 < 23/30-5 < 1-7/30. As a result, combining ([£55]) . 
dHTMD , flSjjgp , (I3T36D . and (f33TD . we obtain 

u;(x) — w(l) 



1 — x 



< ^2/5-5 + £ l/30-5<5/2 + =.l/5-<$ 4. ^7/30 _,_ ^1/30 



L 2 (-l,l) 



< 1/30-5(5/2 



+ e' /lJV + E 1 



which provides the required result for 5 < 1/45. Combining all together, we proved that C(e) — > 
as e — > in bound (|3.15p . According to the previous construction, this finishes the proof of 
Theorem I 
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3.4 Convergence rate of eigenvalues of A e 

To prove the convergence rate of the Main Theorem, we write the eigenvalue problem A e w = fxuu 
as the generalized eigenvalue problem 

L e _w = 7£~ 2 (L^rV (3.38) 

where 7 = Let us first introduce some notations. For any integer n ^ 1, let w n be an 

eigenvector of Aq for the eigenvalue [i n = 2n (n+i) > anc ^ ^ Un = 2(i-x' i ) ' According to the results 
of section T3.21 w n is identically equal to outside of the interval (—1,1) and its restriction to 
(— 1, 1) is a polynomial which vanishes at the endpoints ±1. In particular, u n S L 2 {M). Moreover, 
u n solves the equation 

-O x +po) U n = H n U n , 



2(1 - x 2 



which means that fi n is an eigenvalue of Aq , with associated eigenvector u n . Conversely, if u G 1? 
is an eigenvector of ^4q for an eigenvalue fi, then w = 2(1 — x 2 )u defines an eigenvector of Aq for 
the same eigenvalue fi. Therefore Aq and Aq have the same eigenvalues {fJ> n }n^i- Similarly, for 
e > 0, A £ and A* have the same eigenvalues {fJ>n,e}n^i, and w nt£ £ L 2 is an eigenvector of A e for 
an eigenvalue [i n , e if and only if u nfi = L e _w n ^ £ is an eigenvector of A* for the same eigenvalue 
/j, n ,e- For convenience, w n and u n are normalized by 

IWIl2( R ) = 1. 

Then, according to Remark 13.71 for any n ^ 1 and any e > 0, we can define an eigenvector u„ i£ 
of A* for the eigenvalue /%, e , in such a way that 

w n ,£ — ► M n in L 2 (M) as e — > 0. 

We also define 

w n,e := A t Tl)£ (^l) 1*n,e = £ ^%U n ^ e . 

Then, we have the following lemma, which gives directly the rate of convergence of 7 nj£ = l/fx ne 
to 7 n = 1//U n hi the Main Theorem. 

Lemma 3.9 Let m, n ^ 1 6e too integers and fix 5 > small. The following alternative is true: 

• Ifm^n, then | J_ x w n n mi£ dx| < e 1/3_<5 . 

• Ifm = n, then \ j ] _ l w n u mt£ dx\ > 1 and - fj, n \ < e 1/3 ~ 6 . 

Proof. We prefer to work with 7„ )£ = l/fi n ,e and 7 n = l/fj, n - The eigenvector of ^4 e , w m ,£ = 
'yf n A e w rrije solves the problem 

- w m,e( X ) = lm u m,e, -1 < X < 1, 

while the eigenvector w n = ^uAqWu solves the second-order differential equation 

-2(1 - X 2 )w'n(x) = J n W n (x), -1 < X < 1. 
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Multiplying the first equation by w n and integrating by parts on [— 1 + e 2 ^ 3 , 1 — e 2//3 ], we obtain 



(7m - In) / 'i',, 'hn 

l\x\<l-e 2 / 3 

where 



>£ dx = [w' n w mt£ - w n w' m ^ | 



'"■-J x-=-l+e 2 / 3 ^ n / w n@m,edx, (3.39) 



Wm,e{x) 



2(1 -X 2 ) 

By Lemma [2T2l since H-Liu^^H^ = 7mll n m,e||L 2 — > 7 m as e — > 0, we obtain 

ll w m,ellL°°(l-e 2 / 3 <|x|<l) ^ H w m,elU 00 (K) ~ (3.40) 

\\ w me\\L°°(i-t*/*<\x\<l) < \w mt£ (-l)\ + \w mt£ (l)\ +£ 2/3 ||^ ie || LO o (1 _ e 2/3 <N<1) < e 2/3 .(3.41) 
The last term in the right-hand-side of (|3.39p is estimated by 



w n e m;£ dx 



II^WIlL^lzKl-e 2 / 3 )- 



'H<l-£ 2 / 3 

The function 6 m>£ (x) solves the second-order differential equation for |x| < 1 — e 2 / 3 : 
2a» r„\ i Jl\a r„\ _ ^ n r„\ T „u_ „ („\ _ w m,e{ x ) 



(3.42) 



— e 



E (x) + 2(1 - x )0 m , £ (x) = e g m>£ (x), where g m>£ (x) 



2(1 - x 



2\ ' 



We infer that 

| 5m , e (±(l - £ 2 / 3 ))| < 1, |<4, £ (±(1 - , 2 / 3 ))| < e- 2 '\ 
We take a scalar product of (|3.43p with Q m ,e and obtain the bound 

e2 ||0m,ellL2(|a;|<;l- e 2/3) + e2 ^ 3 |l^m,e||^2(| a .| <1 _ E 2/3) % £ 2 |^m,e(l ~ £ 2/3 )l \0' m ,e( l ~ 

+£ 2 \0m, £ (-l + e 2/3 )\\9' m ^ e (-l + £ 2/3 )| + £ 2 ||6'm,e|lL2(| I .| <1 _ e 2/3)||5m,ellL 2 (|x|<l-e 2 / 3 )- 

By Lemma 12.51 for a = 2/3, we have for any small 5 > 

\u nh£ {±(l-e^))\= £ -\(L%)- l Wmt£ )(±{l-e 2 l 3 ))\ < s~ 5 , 



(3.43) 



(3.44) 



ii 



'(±(l-e^))\=e-\(L%)^w m J(±(l-e 2 / 3 ))\ < . 



.-8 



The bounds (pOIjh (pOHjh and (13^471) . induce, if J < 1, 

|0 ro , e (±(l - e 2/3 ))| < \u m , £ (±(l - e 2 / 3 ))| + bm, £ (±(l - e 2/3 ))| < 

IC, e (±(i - e 2/3 ))l < l< e (±(i - e 2/3 ))l + W m ,eW ~ e2/3 ))l £ £_2/3 " 5 - 

On the other hand, it follows from the definition of g m ,e m 

(g35D that for x € (-l + e 2/3 , 1 
<, £ W = 2(1 - x 2 )g'^ e (x) - &xg' m £ {x) - 4g m , £ (x). 



(3.45) 



(3.46) 
(3.47) 



(3.48) 
(3.49) 

-e 2/3 ), 
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We multiply this identity by g'^ and integrate over (-l + e 2 / 3 ,l-e 2 / 3 ). We get 

-is 2 ' 3 



1-e 2 / 3 

(l-X 2 )\g'^dx + 8 

-1+e 3 / 3 

1-e 2 / 3 



-1+e 2 / 3 



^ 2/3 

Wm, e g'L,e dx + 4 [ x 9m,e( x ) 2 + 9m,e( x ) 9m,e( X )) _^+ £ 2/3 ) 



l+£ 2 / 3 



which implies thanks to Lemma 12. 11 (|3.44|) and the Cauchy-Schwarz inequality 



,£»L 2 {-l+e 2 / 3 ,l-e 2 / 3 ) 



+ IIS'm,ell_L2(_ 1+e 2/3 jl _ e 2/3) ^ \\9rn,s\\L 2 (-l+e 2 / 3 ,l~e 2 / 3 ) + £ 



It follows that there exists C > such that 



e2/3 (ll5m,elL 2 (-l+£ 2 / 3 ,l-£ 2 / 3 ) _ Ce 2/3 ) + ll5m,ellL 2 (~l+e 2 / 3 ,l-e 2 / 3 ) 

As a result, 



llflm,elli, 2 (|a;|<l-e 2 / 3 ) ~ £ ll#m,e lll, 2 (|a;|<l--e 2 / 3 ) 



< 



4/3 .(3.50) 
< e- 4 '\ (3.51) 
(3.52) 



Then, thanks to (pT4"5]) . ([333), ([S32D and (13T521 . we obtain 

12 < 



£ ll^m,ellz, 2 (|a;|<l-e 2 / 3 ) 



3-, + £ 2/3 ||6' mi£ || 2 2( | x | <1 _ e2/ 3 ) < e||0m,e|li2(|a!|<i_ e 3/3) 



+ £ 4/3-25_ 



Therefore, there exists e-independent constant C > such that 



£2 H^m,ellz, 2 (|j;|<l-£ 2 / 3 ) 



+ £2/3 (ll0 m , 



Il 2 (|:t|<1-£ 2 / 3 ) 



C£ l/3\ < £ 4/3- 



2<5 



Thus, 



ll 6 'm,ellL 2 (|x|<l-e 2 / 3 ) ~ £ & 



We deduce from dCTjl . d320|) . d33I|) . d322|) and ([3331) that 

|-l-£ 2 / 3 
/_ 1+£ 2/3 



(7m - 7n) 



< =-1/3-5 



(3.53) 



(3.54) 



,l_ £ 2/3 



If m / n, then - 7 n | > 1 and therefore J_ 1 + e 2/ 3 w n u m<£ dx 
L 2 (M), using the Cauchy-Schwarz inequality, we obtain 



< e 1/3 5 . Since u m , £ -> u m in 







^l-£ 2 / 3 




y w n u my£ dx 




/ w n u m:£ dx 


+ 






J-1+E 2 / 3 





l-e 2 /3<[ a; [ < i 

which is the estimate of the first alternative. If m = n, since u. 
l[_ 1+£ 2/3 jl _ e 2/3]n nj£ -> u n in L 2 (M), and thus 



< e l/3-^ +e l/3|| ? 



n.e 7 "ti 



in L 2 (R), we also have 



1-e 2 / 3 

w n u n , £ dx — > 

, 1+£ 2/3 e^O 



w n u n dx 



-1 



1 w 2 , 



! 2(1 -X 2 



-(fa > 0. 



Combined with (|3.54[) . it gives \ j n ,e — 7n| ^5 5 > which is the second alternative. 
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4 Eigenvalues of the spectral problem (11.51) 



As we have seen before, if (u, w) G L 2 (IR) x L 2 (M) solves system (ll.5p . then u> is an eigenvector 
of A £ associated to the eigenvalue I/7, where 7 = — A 2 /e 2 . In other words, to solves the two 
fourth-order differential equations 

e 2 (-d 2 + ^(x 2 - l)) 2 w(x) = jw(x), for \x\ > 1, , 4 ^ 

-2(1 - x 2 )u/'(x) +£V"(x) =710(0;), for |x| < 1, l ' J 



which also means that u; solves the generalized eigenvalue problem (|3.38p . Since w E L 
we have (L^) _1 u> G flj 2 c (R) C C X (IR) for any fixed e > 0. Prom the generalized eigenvalue 
problem (I3.38p . we infer that w is twice continuously differentiable on R and w"'(x) has jump 
discontinuities at x = ±1: 



i//\x=l+o 2 /i\ ni\x=— (l— 0) ^ 

w L=i-o = ^(i). ™ L=-(i+oj = ( 4 - 2 ) 



Solutions of the first equation of system (|4.ip on the outer intervals {|x| > 1} can be constructed 
analytically. Solutions of the second equation of system (|4.ip on the inner interval (—1, 1) can be 
approximated numerically. Following to a classical shooting method, we shall find numerically 
an estimate on the convergence rate of jn^ to 7 n as e —* 0, for a fixed n ^ 1. The convergence 
rate we observe numerically is faster that the one in the Main Theorem. 

For convenience, we will only consider even eigenfunctions w(x) near 72m— 1 = 4m(2m — 1) 
for an integer m ^ 1. A similar analysis can be developed for odd eigenfunctions near 72™ = 
Am{2m + 1) for an integer m ^ 1. 

4.1 Asymptotic solutions on the outer interval 

For a fixed value of 7 > 0, w solves the first equation of system ()4.ip on [1, +00) if and only if 



£ 2 M * ' e 2 



„ x 2 - (1 -eJj)\ ( x 2 -(l + e % /7)\ 
-dl + \ 2 Vn -a 2 + 1 Vn ) w. (4.3) 



Thus, linear combinations of solutions of the second-order differential equations 

= (-g> + *-v+^ y (4.4) 

for v = ±^/7 provide solutions of the fourth-order differential equation (|4.3p . We shall see 
that they are the only solutions of (|4.3p . First, the following lemma gives a set of two linearly 
independent solutions of (|4.4j) . 



Lemma 4.1 Fix !/ £ I. There exists a constant C > suc/i i/iai /or e > sufficiently small, 
the equation 

-tf'(x) + 7 ^ V(x) = -V(x), x^l (4.5) 
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has two linearly independent solutions ijff and ipg £ such that for x ^ 

r/(VTT^(l + *)) = a(x)Ai f il±^!iM j (i + Q^(*))) , 

^(VT+^(i + x)) = a(x)m ( ^ 1 + £ ^ {x) ^J (1 + Q^{x)) , 

u>/iere £(x) := f| J7; -y/i(2 + t)dt) ^ , a{x) := (^'(x)) -1 / 2 and Q^ £ , Q# e satisfy the bound 



I|QTIU°°(K + ) + IIQB e |U°°(iR+) ^ Ce 2 / 3 . 

Moreover, 

(VgW = 2 1 / 3 Ai-( £ V3 2 - 2 /3, ) _ 6 i/ 3r(2/3) 

^ £ (!) e 2 / 3 Ai(e 1 /32-2/3 iy ) ^ + u ^ ^" e 2/3 r(1/3 ) ^ + - 

where C(e 1//3 ) and C(e 2//3 ) in J^.6| ) are uniform inv^K, for any compact set K C R. 

Proof. See Appendix IA.31 



(4.6) 



Remark 4.2 ./Vote i/iai solutions of ^4-4\ ) can be expressed in terms of the Whittaker's functions 
of the parabolic cylinder equation. The connection of these functions with Airy functions, similarly 
as in Lemma \4-l\ was studied by Olver FD) / using asymptotic formal methods. 



Corollary 4.3 Let n ^ 1 and w £ E L 2 (M) be an eigenvector of the generalized eigenvalue problem 
A3. 38\) for the eigenvalue 7ne- Then, there exists constants c+ and c_ such that 

w £ (x) = c + ip^^ ,£ (x) + c-^j±^' £ (x), x > 1. (4.7) 

Moreover, 

M1} = -r(i/3)^<(i) 1/3 ) <( = jWVffldS (i + 0(£ v 3) ) , (4 . 8) 

W 6 1 /3r(2/3) V v 7' eU 6 1 /3r(2/3) V v V v 7 

Proof. First, we remark that if 7 > 0, then yy^ ' ' ,tp^ ' and ip B ' are four linearly 
independent solutions of the fourth-order equation (|4.3p . Indeed, if Cj[,Cg are constants such 
that 

<%*f " + C^f £ + C~ A ^ £ + C~ B ^ e = 0, (4.9) 
applying the operator — <9 2 + ^-?r- to (|4.9|) . we obtain 



Combined with (|4.9p . it gives 



C+^' £ + C+^' £ = and CJ^^' £ + C-^ B ^ ,£ 
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Prom Lemma [4. II and from the asymptotic behaviour (|2.38p of Ai and Bi, we deduce that for any 
f£K, tp 1 ^ 6 and ip 1 ^ 6 are linearly independent. As a result, C\ = Cg = = Cg = 0. It follows 
that the only solutions of (|4.3p which vanish at infinity, are the linear combinations of VOi ^ 

and ^ j4 v/ ^' e . It results in the decomposition (|4.7p . Since 7„ je — ► 7 n as e — > 0, the asymptotic 
expansions (|4.8p come from (|4.6|) and the identities 



<(1) 
<(l + 0) 



c + ^f^' £ (l)+c_^^' £ (l), 



^'(Tn,,) 1 / 2 
^^(Tn,,) 1 / 2 



-^^(lj + c^^il) 



-c + (^f^' £ )'(l) + c_(^^ £ )'(l) 



+2e 



-2 



c+V'* (1) + c_V 



u£'(1-0) + 2e" 



(1) 



(1)+C_£ 



(1) 



Remark 4.4 Asymptotic limit (f^.6| ) implies that for < £ < 1, i/te eigenvalue A^j o/ i/ie se//- 
adjoint problem L e _w £ = \ e n w £ satisfies a sharp bound 

C~e 2 / 3 < |A£ - A n | < C+e 2 / 3 (4.10) 

/or a fixed integer n ^ 1, where X n = - 2 > A^ i> s the n th eigenvalue of L e _ and < C~ < < oo 
are some constants. Indeed, differential equation L £ _w = Xw has analytic solutions for even 
eigen] 'unctions 

for \x\ < 1, 
for \x\ > 1, 



w 




where c is a constant. Notice that for A > fixed, v = eX stays in a compact subset of K when 
e goes to 0. Continuity of w(x) and w'(x) across 1 leads to an algebraic system, where c can be 
eliminated and X is found from the transcendental equation 



cos(y / A) _ ^f £ (l) ^ 2/3 _ r(l/3) 



VAsin(^A) " (^ A,e )'(l)e- 6V3r(2/3)' 

where we have used J^.6p . We deduce that for some integer m ^ 1, vA = y / A| m _i = \/X~2m^i — 
6 m (e), where V / A| m _ 1 = 7r ( 2r ^~ 1 ) formal are the roots o/cosV% and8 m {e) ~ £ 2/3 (2m^-i^r^i/3) 
proves ^Tffl for n odd. For odd eigenf unctions (n even), the analysis is similar. 



4.2 Numerical solutions on the inner interval 

Unfortunately, Remark 14. 4 1 is not useful in the context of the non-self-adjoint system 114. IB because 
we do not know explicit analytic solutions of the second equation of system (|4.ip . Therefore, we 
use a numerical method to approximate these solutions on the inner interval [—1,1]. 
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Considering even eigenfunctions of (|3,38p we let W\{x) and W2(x) be two particular solutions 
of the second equation in (|4.1|) on [0, 1] subject to the boundary conditions 

" ioi(l) = 1, w'((l) = 0, 1^(0) = 0, <(0) = 0, 
«*(!) = 0, w' 2 '(l) = l, 4(0) = 0, <(0) = 0. 



Then, a general even solution of the second equation of system (|4.ip writes 

aiwi(x) + a,2W2{x), < x < 1, 



w[x 



(4.11) 



for some constants a%, 0,2- The continuity of w{x) and w"(x) across x = 1 leads to the scattering 
map from (01,02) to (c+,c_) in the solutions (|4.7p and (|4.1ip . which is solved uniquely by 



ai =F £7 x / 2 a2 



2^ 



where for conciseness, 7 nj£ is simply denoted 7. The continuity of w'(x) and the jump condition 
(14, 2p on w"'(x) across x = 1 lead to a linear system on (01,02) in the form 



r 2/3 



tfli(l) 



Ol + 



£7 



-1/2 r/ 



r 2/3 



Ol + 



where 



r2/3r,/,^' £ 



WOT 7(1) 



2^ 



+ 



1) 



(1) 



et/„ 



£4 



4(1) 



r 5/3 <(l) 



«2 
«2 



0, 
0. 



2/3 Wf' £ )'(i) , £ 2/3 te^ ,£ )'(i) 



2^ 



+ 



1) 



1) 



■>a w "^a w *ya 2 ^a 

By the ODE theory, unique classical solutions u>i(x) and W2(x) exist for any e > and the 
dependence of uii^x) on e is analytic for e > 0. If there exists a simple root of the determinant 
of the linear system for a particular value £0 > 0, the root persists for other values of e > near 
£ = £o- This method is used for tracing eigenvalues 7(e) of the spectral problem (|3.38p as e — > 0. 

To do it numerically, we approximate solutions w\{x) and 1^2 (#) with the second-order 
central-difference method on a uniform grid with the grid size h = 0.005. The numerical method 
is explained in Appendix IA.51 On the other hand, the values of U p and U m can be evaluated 
from the asymptotic formula (|4.6f) for e £ [10 -6 , 10 -4 ] with 20 data points. Using these approx- 
imations, the determinant of the linear system for (01,02) is plotted versus 7 near 7 = 71 = 4 
and 7 = 73 = 24 and its zero is detected numerically. Then, the zero is plotted versus e and its 
best power fit is used to detect the convergence rate of I7 — 7 n | ~ Ce v . The numerical zeros and 
the best power fit is shown on Figure 3 for 71 = 4 (left) and 73 = 24 (right), while the numerical 
approximations of the eigenfunctions for e = 10 -4 are shown on Figure 4 (dots) together with the 



limiting profiles obtained from the polynomial C 2 (x) and C 4 ±/z (x) at e = (dashed lines). 
The numerical values of the power of the best power fit are found to be 1.9959 for 71 
1.9662 for 73 = 24, which suggests that the sharp asymptotic bound is 



,-1/2, 



4 and 



l7n,< 



7n 



< 



for n 1. Finally, Figure 5 shows the ratio 01/02 obtained from the linear system for e = 10 -6 
in 7 near 71 = 4 (left) and the values of the ratio at the non-zero solution of the linear system in 
£ (right). The power fit was found to be 1.99998 and it illustrates that lim e _^o a i/ a 2 = 0, such 
that lim e ^o w{x) = W2(x) (up to renormalization) . 
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Figure 3: The numerical zeros of the determinant of the linear system (dots) 
and its best power fit (dashed line) for 71 = 4 (left) and 73 = 24 (right). 




Figure 4: The numerical approximation of even eigenfunctions (dots) for 
e = 10 -4 near 71 = 4 (left) and 73 = 24 (right) and the even polynomial 
solutions for e = (dashed line). 




Figure 5: The ratio a\/a2 for the two equations in the linear system versus 
7 for £ = 10™ 6 near 7 = 71 = 4 (left) and for the solution of the linear 
system versus e (right). The best power fit is shown by dashed line. 
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A Appendix 

A.l Proof of Lemma 12.11 

Let us denote by Ai(LL) the smallest eigenvalue of L e _. We first show that Ai (LI) > 1. Let 
X G C C °°(IR) be such that < \ < 1, supp(x) C (-3,3), and \ = 1 on (-2,2). Let 5 > to be 
fixed later (independently of e) . The Max-Min principle ensures that 



vED(L e _) \\V\ 



L 2 



inf \\\v'\\ 2 ,2+ [ p £ \v\ 2 dx\ = min{A (1) , A (2) }, (A.l) 

veQ(L±),\\v\\ L 2=l I J\*\>1 / 



where 



A (1) = inf \\v'\\h + / Pe\v\ 2 dx 

v e Q(LL), \\v\\tf = 1, \ y| a |>i . 

/ w>2 M a <fc>a 

A (2) = inf [ ||i/||^a + / p £ |v| 2 dx 

« £ QC-^Di \\v\\ L 2 = i, V 7|x|>i 



If ||w||i2 = 1 and Ji x i >2 \v\ 2 dx ^ 5, then 

I\ x \>2^ 2 ~ l Mdx > *!\ x>2 \v\ 2 dx > 35. 

Therefore for e ^ 1, 



A« > ^ > 3& (A.2) 



On the other side, let us now take v G Q{L e _) such that ||u||x,2 = 1 and Jj x | >2 M 2 ^x ^ 5. Then 

/ (x 2 - l)\ X v\ 2 dx < / (x 2 - l)H 2 dx, (A.3) 



'H>1 J|s|>l 
and since x'C^) is supported in {2 ^ |x| ^ 3}, we also have in this case 



\(xv)'\ 2 dx 



X 2 \v'\ 2 + 2 X x'vv' + X ' 2 \v\ 2 



dx 

< ll^llz,2( R ) + 2 \\ v '\\l2(R)\\X \\l°°(R)\\v\\ L 2 {\x\>2) + llx'llz^iylMlz^i^a) 

< 2||t;'|| 2 2(R) + 25|| x '|| 2 0O{R) . (A.l) 

Next, since x = 1 on {|x| ^ 2}, 

\xv\ 2 dx ^ y |t;| 2 dx^l-5. (A.5) 
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Thanks to (|A,3|) . (|A.4|) and (|A.5|) . it turns out that 

>1 I X^l ^ 



. y^g) + 2 ^IIX'H|~(R) + !\ X \>lPe\v?dx 



As a result, using (|A.6p . since (x^)|(— 3,3) G Hq{— 3, 3) for u G ff 1 (IR), 

2 A(2) > ^Hx'Hl^R) | . nf /- 3 l w/ | 2da; + ilx|>iPglH 2 ^ 
1 - & 1 - ^ iue^(-3,3) |u>| 2 ofe 

2<5|| X 'II? - IU./112 



1-0 «i6fl-J(-3,3) 1 1 It? 1 1 2,2 



Thanks to the Poincare inequality, we can now choose 5 G (0, 1) sufficiently small such that 
R s > 0. Then, according to (TO]) . ([Oil and ([ATfjl . 

Ax(XI) ^ min^, (1 ~ 2 ^ V (A.8) 

which provides the estimate Ai(Lf_) > 1 for < e ^ 1. The other estimate \\(L £ _) < 1 is a 
direct consequence of (|A.1|1 and of the Poincare inequality. Indeed, the right hand side in ()A.1|) 
is bounded from above by the infimum of the same quantity, taken over v G L 2 (M) such that 

G #o( -1 ' *) and u |{M>i} = 0- ■ 

A. 2 Proof of Lemma 12.31 

To prove Lemma 12.31 we use the following lemma. 
Lemma A.l For e > 0, 

defines a self-adjoint operator on L 2 (R). The spectrum of L £ is made of a sequence of strictly 
positive eigenvalues increasing to infinity, and the smallest eigenvalue satisfies 



Ai(L e ) «e" 



-4/3 



Proof. The first assertion is straightforward. Thanks to the Max-Min principle, \\{L £ ) is given 
by 



Ai(L E ) = inf (Vll^ + 4 / |x|v 2 ^ , 



\H\l2 = 1 

where 



Q{L £ ) = {vE ff x (R) : |x| 1/2 t> G L 2 
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is the form domain of L £ . If v € L 2 (R) and ||u||x,2 = 1, v can be rewritten as v(x) = hw{h 2 x), 
with h > and to G Q(L £ ), with ||«;||x,2 = 1 and ||k/||z,2 = 1. Moreover, h and w are uniquely 
defined this way, and we have 

n i\ 

II 2 



\v'\\ 2 T ,=h 4 



and 



Thus, 



where 



/ \x\v 2 dx = h 2 I \x\w 2 dx. 
Jr Jr 

X^L') = inf (a 4 + e- 2 ^" 2 /?) = ( JL + 2 1 / 3 ) /^V 4 / 3 , 

/3 := inf / \x\w 2 dx. 

w e Q(L e ) Jr 



iu e Q{L e ) 

\\W\\ L 2 =1,\\W\\ L 2 =1 

The lemma follows if we prove that /3 > 0. Let us assume by contradiction that /? = 0. Let 
( w s)s>o De a minimizing sequence, that is Hit^H^ = ||w^||i,2 = 1 and J R \x\w 2 dx —> as 5 —> 0. 
Let x G be such that < x < 1, supp(x) C [-1,1], and x = 1 on [-1/2,1/2]. For 

a > 0, we also define Xa(a^) = x( x / a )i as well as ws, a '■= Xa.ws- Thanks to the Poincare inequality, 
a : = inf > 0, and then inf -f^ 2 - = 2 > 0. Thus, 

o 2 

lK,Jl£a(R) > -2\\ w S,a\\h(R) 



a* 
a 



2 



a2 IKII L 2 ( _ fif) 



2 IL„ ||2 



» i{ i -lL yA ^ dx )- (a - 9) 

On the other side, since x'( x ) is supported in {1 ^ \x\ ^ l}, we have 

H,a\\h = I ((x'a) 2 W 2 s + 2XaX'aWsw' 5 + X 2 a(w' 5 ) 2 )dx (A.10) 
JR 

< — ^^2°° (E) INHi 2 (|<M<a) + ^IIX / ||L-(R)ll^llL2(f<|x|<a)ll^llL2(R) + 1 1 ^5 1 1 |s (M) - 

According to the assumption, given a > 0, we can find 5(a) sufficiently small such that 



/ 

JR 



x \ w 8(a)d x ^ fl2 - 



Then, 

/ w 2 dx ^ / u> 2 dx ^ - / |a;|tt> 2 dx ^ 2a. (A. 11) 

J%<\x\«l J\x\>^ a JR 
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It follows from (|AJ|, ([AlU]) and (|ATL]) with 5 = 5(a) that 

2||x'lli- fln 2 3 / 2 || X '|| LO 



^(l-2a)^-" /llio 



+ 



+ 1. 



a a 1 / 2 

Letting a go to yields to a contradiction, which completes the proof of the lemma. 



Thanks to the Max-Min principle, we know that the lowest eigenvalue of L £ + is given by 
where 



(A.12) 



L 2 



Q(L e + ) = {v£ H^R) : E L 2 ( 

is the form domain of L e + . The statement of Lemma 12.31 is equivalent to Ai(Lfj 
first prove the upper bound on Ai(L+). Let us define v £ on R as 

x - 1 + e 2 / 3 for 1 - e 2 / 3 < x < 1, 

f e (x) := { — (x - 1 - e 2 / 3 ) for 1 < x < 1 + e 2 / 3 , 
elsewhere. 



and denote q(x) := e 2 q £ (x) = 2(1 — x^l/u^n + (x 2 — l)l/| a .i >1 ').. Then 

ll u elll 2 (R) = 2e2 ^ 3 ) IK 
and since q(x) ^ 4|x — 1| for |x — 1| ^ 1, 



2f_ 
3 



q £ \v E \ dx ^ _ 2 



l+e 2 / 3 



Jl_ e 2/3 



|1 — x\v e dx 



2£ 2/3 



As a result, 



v +y 2e 2 /3 



It remains to find a bound on Ai(Ll) from below. Let us first introduce the two intervals 



-D+ := < x > 0, ofx) ^ 



1 



v 7 ^ /3 
^'V2 



, D_ := 1x^0, q(x) < 



and denote D := D+ UD-. If v G Q(L+), [|«[|l 2 = 1 an d Jd ^ 1 — e 1 ^ 2 , then 

1 



v\ dx ^ / dx ^ 



£ l/2 

M 2 cix > > 4e 2 / 3 



for sufficiently small e > 0. As a result, thanks to (|A.12p and the upper bound on Ai(L^), we 
deduce that 



Ai(£ e H 



inf 

v 6 Q(L%), 

\M L 2=1, 

J D \v\ 2 dx ^ 1 - e 1 ' 2 



J 



+ / Qe\v\ 2 dx 



(A.13) 
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Prom now on, we assume that v £ Q(L e + ), \\v\\ L 2 = 1 and J D \v\ 2 dx ^ 1 — e 1 ! 2 . Let x £ 

be such that ^ x < 1, supp(x) C [-1/2, 1/2] C M\D, and x(x) = 1 for z G [-1/4, 1/4]. We 

also define p := 1 — x- In particular, p = 1 on D, thus 



2 > 

L 2 ^ 



/ |v| 2 dx > 1 -e 1/2 , / g|H 2 «&< / <?H 2 dx, (A.14) 
J d Jm. Jr 

and since p' is supported in M\D, for some C > 0, we have 



\{pv)'\ 2 dx < IIp'III-^IIwII^^) + llt/Hiafu) + 2||p|| l ^ {r) ||/9 / || l ^(m)||w / || L 2 (r) ||v| 

+ ||«'|||2( R ) + Ce 1/4 ||v'|| L 2 {R) 
< 2(||^|| 2 2(R) +C £ 1 /2 ) . (A . 15) 

Therefore, combining (|A.14|) and (|A.15p . we obtain, for e sufficiently small, 

\\(pvy\\ 2 L i+J R qe\pv\ 2 dx 2(\\v'\\ 2 L2{R) +Ce 1 / 2 ) + f R q £ \v\ 2 dx 



h * 1-^ 1/2 



< 2(\\v'\\ 2 L2 + f q E \v\ 2 dx) + 2Ce 1/2 . (A.16) 
Taking the infimum in v in (jA~16|) . we infer thanks to (fA~T3i) that 

2A l( ^) + 2^/ 2 > inf IIMIIj. + J^H 2 ^ (A1?) 

=1, 

/ D |-u| 2 dx ^ 1- £ V2 
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Therefore, since q(x) ^ 2\x — 1| for x ^ and q(x) ^ 2\x + 1| for x ^ 0, and decomposing 
/to = v i + 1*2 with v\ supported in (— oo, —1/4] and v 2 supported in [1/4, +oo), we have 

2Xi(L E + ) + 2CS 1 ' 2 

> inf \Wi\\h + LQe\vi\ 2 dx+ \\v' 2 \\ 2 L2 + J R q £ \v 2 \ 2 dx 
vi,v 2 £Q(L%), \\vi\\ 2 L 2 + \\v 2 \\ 2 L 2 

supp(ui) C (— oo,— 1/4], 
supp(t> 2 ) C [1/4, +oo) 

> inf HWh + g J R |a + l|M 2 £fc + + | / R |x - 1|H 2 ^ 
»i,"2eQ(£ £ + ), INI| 2 + ||«2|||2 

supp(^i) C (—00,— 1/4], 
supp(ti 2 ) C [1/4, +00) 

> inf HWh + jjg> + i|bi| 2 rfx + |K||| 2 + jjg> - i||u2| 2 rfg 

»i,»2eO(i £ + ) INl£a + IM^a 

HJlg + Jr J R |x|H 2 dx + + ^ / R |x|M 2 <fa; 

»i,"2eg(iy Ikill^ + [|u 2 [| 2 2 

inf Ihjjljg + jr/ R MM 2 cfc + IKII|2 + ^J R |a|M 2 <*c 
ui,« 2 eQ(iy ||^2 + \\V2 ||| a 

> inf / Ibillja + Jr J R Nh| 2 dg + |K|| 2 2 + Jr J R \x\\v 2 \ 2 dx \ 
vi,v 2 e Q(L E + ) y 2 2 J 

\\vi\\ L 2 5? [ML 2 = 1 

> \ ^f f|Klli 2 + 4 / \x\\v 2 \ 2 dx) >l\i(L*)> e-*/ 3 , (A.18) 
2 v 2 e Q(Z^) V e JR J 2 

\\V 2 \\ L 2 = 1 

where we have used Lemma lA.ll in the last estimation. ■ 
A. 3 Proofs of Lemmas 12.61 and 14.11 

Proof of Lemma 12.61 The proof of Lemma 12.61 relies on WKB approximation techniques, 
explained for instance in [Mj . If we define w(x) := — x), it is equivalent for ip to solve (12. 37|) 
or for w to solve 

e 2 w" - 2x(2 - x)w = 0, x € ( 0, - J . (A.19) 



for := ^4 to solve 



2/3 



In the new variable £ = £(x) := ( | y/2t(2 — t)dtj , it is equivalent for w to solve (1A.19P or 



o ef = e 2 *(o«, ee(o,eo), (A.20) 
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where £ := £(3/2), a(x) := (£'(x)) 1/2 , and £(£) := -a"(x)a 3 (x). Next, we look for v in the 
form v(£) = Ai (zi/^j (1 + Q(Q)- Using that Ai(£/e 2 / 3 ) solves the homogeneous equation 

it is equivalent for v to solve (jA.20j) or for Q to solve 

2 



Ai (^) 

By integration, (|A.21[) is equivalent to the integral equation 



ML)A] \4j3) ( 1 + Q(0), £g(0,&). (A.21) 



r& ,r,Ai'-i 



Q(0 = WXO == / ° f ) e2/3 ( 2 d^(77)(l + Qfa))^, (A.22) 



where F maps C°([0, £0]) into itself. A change of variable provides 



F(Q)(£) = e 2 ' 3 (j^ Ai{u) -2 duAi (J?,) 2 ) S (r,)(l + Q( v ))d v . 



Thanks to the asymptotic behavior (|2.38j) . f(x) := /J : Ai(y) 2 dyK\(x) 2 ~ 2y^x ^ — * -\-oo. In 
particular, / is bounded on E + . We deduce that for any £ G (0,£o), 

|(F(Q))(£)| ^£ 2/3 ||/IIl-( K+) f° \Kn)\dr}{l + ||QIU»(o,&))- 
Since 5 is clearly continuous on (0, £0] an d 

9 • 2 2 / 3 

<y(£(z)) — ► 56Q as x -» 0, 
we deduce <5 G L x (0,£o)- Thus, if Q G C°([0,£ ]), then 

||F(Q)|| ioo((Uo) < £ 2/3 ||/||L-(M + )||5|U 1( o lf o)(l + llOIU» ( o, &) ). (A.23) 
Moreover, if Qi,Q2 G C°([0, £0]), we get similarly 

ll^(Qi) - ^(Q 2 )IU~((Uo) < ^ 2/3 ll/llL-(M + )PllLUo, f o)HQi " Q2|U<-(0,&). (A-24) 

From (|A.23|) and (|A.24|) we infer that, if we take C := 2||/||£,oom + )||£||i,i(o ) fo)> f° r e sufficiently 
small (namely e 2 / 3 < 1/2C), F maps the ball of radius Ce 2 ! 3 in C°([0,£o]) into itself, and is a 
contraction on that ball. Then, F has a unique fixed point Q such that ||Q||l°°(o,£o) ^ Ce 2 /' 3 . Such 
a fixed point of F gives a C 2 solution of flASH) on (0,£ ). Defining as Q^(x) := Q(£(l - x)) 
and applying the sequence of substitutions backwards, we found a solution i\) £ A of the system 
(I2.37P with the required bounds. 
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For the existence of the solution tpg, we proceed similarly. Namely, we look for a solution to 
([A720]) in the form v(g) = Bi (-4m) (1 + Q(0)- lt is equivalent for v to solve (|X20|) or for Q to 



solve 



Bi 



r-2/3 



«5(£)Bi 



r-2/3 



(i + Q(0)» ee(O^o)- 



(A.25) 



Since g(x) := Bi(x) 2 J x +0 ° Bi(u) 2 tin ~ ^= as x — ► +00 thanks to the asymptotic behavior (|2.38p 
again, (7 is bounded on K + . It enables us to prove the existence of a fixed point to the functional 
G:C ([0,£o])^C ([(Uo]) defined by 



G{Q)(t) 



i /•£ Bi 



.'J 



J V B: 



( £ 2/3) 



r d*<y(r7)(l + Q(r7))dr;, 



similarly to what has been done for F. 

The linear independence of ip £ A and ip £ B follows from the linear independence of functions Ai 
and Bi. ■ 



Proof of Lemma 14. 1L The proof is very similar to that of Lemma 12.61 so that we will 
only point out the differences. It is equivalent for tp to solve f|4.5j) on + ev, +00) or for 
w(x) := ip(Vl + + x)) to solve 

e 2 w"{x) - x{x + 2)w(x) = (A.26) 

on M + , where e := e/y/1 + ev. We look for w in the form w{x) = a(x)v(£(x)), where = 

(§ £ yjt{2 + t)dt) and a(x) = ^'(x))" 1 / 2 . Then, it is equivalent for w to solve ([AT26]) on M+ 
or for v to solve 

e 2 v"(0-^(0 = e 2 5(Ov(0 (A.27) 

on M + , where the function £ is defined by <5(£(x)) = — a" (x)a(x) 3 . Since a G C°°([0, +00)) 

and <5(f) ~ 7£" 2 /1024, we deduce that 5 G L 1 (M+). Then, the existence of Q G C°(M + ) with 

l|QIU°°(K+) ^ £ 2/3 , such that = Ai(£/e 2 / 3 )(l + Q(£)) solves (|AT27|) . is established like in the 
proof of Lemma [2. 61 applying the fixed point theorem to the functional F defined in (|A.22p . with 
£0 = +00. Therefore, we obtain i^jC- ^he expression for ip 1 ^ is obtained similarly as in Lemma 
12. 6[ Next, the expression of ip^^x) at x = yjl + ev yields 

^/(VTT^) = o(0)Ai(0)(l + Q^(0)) = a(0)Ai(0)(l + 0(e 2/3 )), (A.28) 

and similarly 

= a'(0)Ai(0)(l + 0(e 2/3 )) + a(0)£'(0)Ai'(0) £ - 2 / 3 (l + 0(e 2 / 3 )) + a(0)Ai(0)e'(0)(Q^)'(0) 
= a(0)C'(0)Ai'(0)£- 2 / 3 (l + O(e 2 / 3 )), (A.29) 
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where we have used that 



(QT)'(o) 



Ai(0)- 2 / M(ri/& 3 ) 2 6(r,)(l + Q v £(v))dr, 



o 



^||5|| Ll(M+) (l + 0(e 2 / 3 ))<l. 



At this point, the function i^f has been defined on the interval [y/l + ev, +00). In the case 
v > 0, we extend into a solution of (|4.5p on the interval [1, +00), thanks to the Cauchy-Lipshitz 
Theorem. We denote I u = [yl + sis, 1] if v < 0, I u = [1, vT+~ez7] if > 0. Then, for any sign of 
v, we have 

i^w-^vr+^i s AWfyu-a,,) < e |(V'r) / (v / i+^)i+£ 2 ii(^Tii^(/,) 

< e^ + ellVriU-C/,) ( A - 3 °) 

and, thanks to (|A.30p 



< \r/{ViT^)\+e\\{^)'\\ Loa{Iv) 

< 1 + e||^ e IU°°(^)' 

thus 

U V a £ Wl~ { i„) < I- (A.31) 
From (LA301) . ([Oil and CA28]) it follows that 

^ £ (1) = o(0)Ai(0)(l + 0{e 1 ' 3 ). (A.32) 

Similarly, 



\(r/ni)-(r/)'(VTT^)\ < e\\(r/)lL~(i„)<\\r/\\ L ~(i„)<i, 

and therefore thanks to (1A.29I) . we get 

($J e )'(l) = a(0)e'(0)Ai'(0)e- 2/3 (l + 0(e 2/3 ). (A.33) 
The limit flM]) follows from fQ2l and ([Q3l . since £'(0) = 2 1 / 3 , and because 

Ai(0) = ... \ . , , Ai'(0) 



3 2 /3r(2/3)' v ; 3V3r(l/3)' 

Notice that all the estimates we made in this proof are uniform in v G K, for any fixed compact 
subset K C R. ■ 

A. 4 Proof of Lemma 13.81 

If / G A' and 92 G D(Lx), we have 

I ( L X'f,(p) DiLx y MLx) I < ||/|U'||£x9?llx < ||/|U'|Mb(Ioe)> 

which provides the continuity of Lx- If / G A' and L x >f = 0, then for every (p G D(Lx), 
(f\Lx<p) x > x = 0- We can a PPly this to 93 = L^x, for any x G A and we get that (/, x = 
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for every x £ X. Therefore / = and Lx> is injective. Let us next prove the surjectivity of Lx 

j X 1x )d{L x )',D(Lx) 



Let T £ D(L X )'. / : x (T, L x l x/ n , r ■., n / r \ clearly defines a continuous linear form on X, 



and for every ip £ D(L 



x, 



(Lxif,<P} D (L x y,D(L x ) = (fi L X<f)x>,X = ( T ' L X lL X<p) D{Lx y :D 



(Lx) 



(T,<p) 



D(L x )',D(Lx) ' 



which means that T = Lx'f- Moreover, the application L x ) : D(Lx)' ► X' we have just defined 
is continuous. Indeed, if T G D(Lx)' and x £ X, 



x L x )T,x) x ,^ x 



T, L x x / 



< 
< 



D(L x )',D(Lx) 
1, 



\ T \\D(L x y\\L x x\\d(l x ) 
\ T \\D(L x y(\\x\\x + \\L~ l x\\ D{L) ) 
\ t \\d(l x y\\ x \\x, 



where we have used the continuous embeddings D{L) C X C H, as well as the continuity of 
G Finally, we show that Lx' is an extension of L. Here, we classically identify 

elements of H to elements of X' (resp. D{Lx)') as follows: if / G H, x £ X (resp. T £ H, 

Lp £ D{L X )), (f,x) x >x = (f\ w ) ( res P- ( T ^)d(l x )',d(Lx) = ( T l^))> where ('I') denotes the scalar 
product in H. Thus, if / G D(L) C X C X' , 

(Lx>f,<P) D ( Lx y,D(L x ) = (f' Ll P)x>X = (f\L<P) = ( L fW) = ( L f^}D(L x y,D(Lx) ' 

which means that L X 'f = Lf. ■ 

A. 5 Numerical methods for inner solutions 

We rewrite the fourth-order equation (|4.ip on [0, 1] in the form 

w"(x) = v(x), e 2 v"(x) — 2(1 — x 2 )v(x) = jw(x), < x < 1. 

Using the finite-difference approximation with the second-order central differences [GPJ, the 
system of differential equations is converted into the system of algebraic equations 

Aiw = v, A2V = yw, 

where v, w are n-vectors of v(x),w(x) represented on a discrete grid {xk}^Zl C [0, 1] with xq = 
and xq < x\ < ... < x n _i < x n = 1. Using an equally spaced grid with step size h = 1/n and 
incorporating boundary conditions w'(0) = 0, v'(0) = 0, we obtain n x n matrices A± and A2 in 
the explicit form, where 



-2 


2 


. 


. 





1 


-2 


1 . 


. 








1 


-2 . 


. 
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. 1 
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and Ai = e 2 A\ — 2diag(l — x 2 ). For the first solution wi(x), with w n = 1 and v n = 0, we obtain 
solutions of the finite-difference equations in the form 



1 



(A\ — 7^4 2 ~ 1 ) 




v = 7 A, w, 




where e n is the n th unit vector in M 71 . For the second solution W2(x), with w n = and v n = 1, 
the finite-difference equations are solved in the form 



The values of w'(l) and w"'(l) are obtained from the three-point finite-difference approximations 



which preserves the second-order accuracy of the numerical method [GP]. 
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